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We show that the emergence of the axial anomaly is a universal phenomenon for a generic three 
dimensional metal in the presence of parallel electric (E) and magnetic ( B ) fields. In contrast to the 
expectations of the classical theory of magnetotransport, this intrinsically quantum mechanical phe¬ 
nomenon gives rise to the longitudinal magnetoresistance for any three dimensional metal. However, 
the emergence of the axial anomaly does not guarantee the existence of negative longitudinal magne¬ 
toresistance. We show this through an explicit calculation of the longitudinal magnetoconductivity 
in the quantum limit using the Boltzmann equation, for both short-range neutral and long-range 
ionic impurity scattering processes. We demonstrate that the ionic scattering contributes a large 
positive magnetoconductivity oc B 2 in the quantum limit, which can cause a strong negative mag¬ 
netoresistance for any three dimensional or quasi-two dimensional metal. In contrast, the finite 
range neutral impurities and zero range point impurities can lead to both positive and negative lon¬ 
gitudinal magnetoresistance depending on the underlying band structure. In the presence of both 
neutral and ionic impurities, the longitudinal magnetoresistance of a generic metal in the quantum 
limit initially becomes negative, and ultimately becomes positive after passing through a minimum. 

We discuss in detail the qualitative agreement between our theory and recent observations of neg¬ 
ative longitudinal magnetoresistance in Weyl semimetals TaAs and TaP, Dirac semimetals NasBi, 

Bii-sSb*, and ZrTes, and quasi-two dimensional metals PdCoC> 2 , a-(BEDT-TTF) 2 l 3 which do not 
possess any bulk three dimensional Dirac or Weyl quasiparticles. 


I. INTRODUCTION 

The axial anomaly plays a profound role in the con¬ 
struction of effective relativistic field theories in odd 
spatial dimensions, involving massless Dirac and Weyl 
fermions^. When massless Dirac fermions in one spa¬ 
tial dimension are coupled to an external electric field 
(E), the separate number conservation laws of the right 
and the left handed fermions are broken by a quantum 
mechanical effect—^, which is described by 

d^(jn,R - 3 p,l) = dtfn ,5 = (!) 

where Nf is the number of flavors, e is the electric charge, 
and j^R/L describe the charge and the current opera¬ 
tors of the chiral fermions. Similarly, in three spatial 
dimensions, in the presence of parallel electric (E) and 
magnetic (B) fields, the anomalous violation of the axial 
current conservation law is described by the celebrated 
Adler-Bell-Jackiw (ABJ) anomaly equation^ 

a »^=SS E - B - (2) 

The emergence of the axial anomaly is intimately tied 
to the topological properties of the Dirac operator with 
unbounded linear dispersion. It is related to the spectral 
flow of the right and left handed states in the presence of 
a topologically nontrivial configuration of the gauge field: 
one sinks below zero energy and simultaneously the other 
one rises above zero energy through an infinite Dirac sea. 
For constructing chiral gauge theories such as the elec- 
troweak theory and the standard model, it is important 


to cancel the anomaly contributions of different species 
or flavors of chiral fermions^. However, it is not exactly 
clear how the notion of the axial anomaly can have any 
impact on the physical properties of a condensed mat¬ 
ter system endowed with a bounded dispersion relation, 
since the infinite nature of the Dirac sea is crucial to the 
existence of the ABJ anomaly in quantum field theories. 

Nevertheless in a seminal paper, Nielsen and Ni- 
nomiya— first proposed that the axial anomaly can give 
rise to a very large longitudinal magnetoconductivity 
(LMC) or a negative longitudinal magnetoresistance 
(LMR) for certain gapless semiconductors, which pos¬ 
sess a linear touching of conduction and valence bands 
at isolated points of the Brillouin zone. Depending on 
whether the bands are nondegenerate or two-fold degen¬ 
erate, the gapless semiconductor is respectively known as 
a Weyl semimetal (WSM) or a Dirac semimetal (DSM). 
The recent discovery of a topological DSM—= 2) in 
Cd.aAsaA— • Na 3 B i 11 i 12 , and the search for a WSM in the 
magnetically ordered (all in-all out) 227 pyrochlore iri- 
datesi^, have rekindled the interest of the condensed mat¬ 
ter physics community in experimental verification of the 
axial anomaly. Very recently a WSM phase with an even 
number of left and right handed pairs has been discovered 
in noncentrosymmetric TaAsiA—. The ideas are being 
put forward that the verification of the axial anomaly in¬ 
duced negative LMR can serve as a diagnostic tool for 
the existence of WSM and DSMiA—. A DSM phase 
with an odd number of flavors also arises at the quantum 
phase transition transition between a strong Zi topo¬ 
logical insulator and a trivial band insulator in various 
materials^ such as Bii-^Sba^A— , BiT^Si-aSe ^ 27 ' 28 , 
(Bii_ 3 ,In a ,) 2 Se 3 29 i 30 . Similarly, a DSM with even num- 
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ber of cones are known to appear at the transition be¬ 
tween a crystalline topological insulator and a trivial in¬ 
sulator in Pbi_ x Sn x Te2i~« 3 . A different kind of gapless 
fermion, known as a massless Kane fermion, is realized in 
Hgi_ x Cd x Te, at the band inversion transition between a 
parabolic semimetal with quadratic band touching (Hg 
rich side) and a semiconductor (Cd rich side) 31 ' 34 . The 
gapless Kane fermion is constructed out of three Kramers 
degenerate bands, and at the degeneracy point there is 
also a parabolic hole like band in addition to the four 
component massless Dirac fermion^. 

Recent magnetotransport experiments on Bii_ x Sb x , 
with x tuned at the band inversion transition^, and 
ZrTeg2£ have identified a large negative LMR. At very 
low fields the LMR is positive, which eventually be¬ 
comes negative with increasing magnetic field strength. 
In contrast, when E and B are orthogonal, these mate¬ 
rials show strictly positive transverse magnetoresistance 
(TMR). The positive LMR for weak magnetic fields is 
associated with the weak antilocalization effects due to 
the presence of strong spin orbit coupling in these mate- 
rials^I. The analysis of negative LMR data leads to the 
existence of a large positive LMC varying as B 2 . In these 
experiments, the concomitant observations of negative 
LMR and positive TMR have been claimed to be the ex¬ 
perimental verification of the axial anomaly of a WSM^ 
(for Bii_ x Sb x ) or a DSM2£ (for ZrTes). Independent of 
whether the axial anomaly is causing the negative LMR 
or not, the very observation of a large LMR does not fit 
into the theme of classical theory of magnetotransport. 
A classical theory predicts zero LMR, due to the absence 
of a Lorentz force when E and B are parallel^, and the 
experiments are clearly pointing toward the existence of 
a subtle quantum mechanical effect. It is also important 
to note that the negative LMR is not caused by the field 
induced suppression of the scattering from the magnetic 
impurities. Otherwise, even the TMR will decrease as 
a function of the magnetic field strength, which is not 
observed in these experiments. The existence of nega¬ 
tive LMR, without the notion of axial anomaly was first 
considered by Adam and Argyres for a three dimensional 
electron gas sixty years ago^. 

However, our particular interest in the potential effects 
of the axial anomaly on the LMR has been stimulated by 
the recent magnetotransport measurements on a quasi- 
two dimensional metal PdCoO^- The Fermi surface of 
this material is a corrugated cylinder, composed of a sin¬ 
gle sheet. The ab initio calculations^ and the ARPES 
measurements^ are in good agreement, and do not show 
any evidence for the existence of an underlying three 
dimensional Dirac like band structure. Nevertheless, a 
large negative LMR and a large positive TMR have been 
concomitantly observed, when the field is applied along 
the crystallographic c axis or the interlayer direction. In 
contrast to the two Dirac materials discussed above (i.e., 
Bii_ x Sb x and ZrTes), this delafossite compound demon¬ 
strates a negative LMR for weak magnetic fields. The 
negative LMR becomes more pronounced at low temper¬ 


atures. Its size and the existence over a wide range of 
temperatures (T < 200AT) and magnetic field strengths 
(B < 30T) can not be reconciled with small weak lo¬ 
calization corrections, which usually occur at weak fields 
and low temperatures^. In addition, if weak localization 
is the underlying mechanism, then it will also lead to 
a negative TMR, in contradiction with the experimen¬ 
tal observations. The field dependence of the negative 
LMR over the wide temperature range can be fitted with 
a scaling function of oj c t according to the Kohler rule2&. 
Here, lo c is the cyclotron frequency and r is the tem¬ 
perature dependent transport lifetime in the absence of 
the magnetic field. The data does not suggest a simple 
B 2 behavior of the LMC. These observations have led 
the authors to propose the axial anomaly as the driv¬ 
ing mechanism behind the observed negative LMR. In 
this context, the axial anomaly has been argued to be a 
ubiquitous feature of any three dimensional metal in the 
presence of a quantizing magnetic field, which arises from 
the quasi one dimensional conduction channels caused by 
the partially filled Landau levels (LLs) dispersing along 
the applied magnetic field. This qualitative picture has 
been further supported by the angle dependent magne¬ 
toresistance oscillation (AMRO) measurements. When 
the magnetic field is tilted with respect to the c axis 
at certain special angles, known as Yamaji angle s 44 ' 45 , 
the LLs become non-dispersive and the material displays 
an extremely large and strictly positive interlayer mag¬ 
netoresistance. As the current is measured along the c 
axis, which is now at an angle with B, the measured 
magnetoresistance is a combination of LMR and TMR. 
When the B field is slightly tilted away from the Yamaji 
angle, the LLs reacquire their dispersion and a negative 
component of MR emerges. More strikingly, the nega¬ 
tive component overwhelms the existing large and posi¬ 
tive component with increasing field strength. This gives 
strong support to the idea that the axial anomaly arising 
in the presence of the one-dimensional Fermi points is 
responsible for the observed negative LMR. 

Motivated by these intriguing experimental results, 
here we develop a comprehensive theory of longitudinal 
magnetotransport in the strong field limit, lu c t 1 . Our 
main achievements are as follows: 

1. We first establish how the dispersing LLs for a 
generic metal can lead to the axial anomaly as 
described by Eq. ©. One may be surprised by 
the fact that we are not placing any stringent 
constraint on the underlying band structure, i.e., 
the system does not need to be a WSM or a DSM 
by any means! This universality is brought out by 
the quantizing magnetic field, which is not a weak 
perturbation. Irrespective of the underlying band 
structure, the spectrum is bunched into highly 
degenerate LLs dispersing along B. We show 
how this field induced dimensional reduction (for 
the purpose of longitudinal magnetotransport) or 
quasi-one dimensionality is ultimately responsible 
for the emergence of the axial anomaly. 
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2. Through an explicit calculation we show how 
the axial anomaly as a quantum mechanical 
phenomenon generically leads to LMR, in sharp 
contrast with the predictions of the classical theory 
of magnetotransport. We employ a Boltzmann 
equation approach to calculate the LMC in the 
presence of short range scattering due to neutral 
impurities and point defects as well as long range 
scattering due to ionic impurities. In the quantum 
limit (B > Bo), when only the lowest Landau level 
(LLL) is partially filled, the ionic impurities make 
a positive contribution to the LMC oc B 2 for any 
metal. In contrast, the short range impurities can 
cause some nonmonotonic behavior of the LMC 
up to a threshold B t and an eventual decrease 
of LMC in the asymptotic limit B —> oo for a 
conventional metal. The decrease of LMC for a 
Gaussian impurity and a point defect vary as B~ x 
and B ~ 2 respectively. Therefore, in the presence of 
ionic impurities alone, the axial anomaly leads to a 
large unsaturated negative LMR, and a LMC oc B 2 
in the quantum limit. This invalidates a common 
assumption in the literature that the axial anomaly 
always leads to a B-linear positive LMC in the 
quantum limit j 18 i 19 i 21 i 35 i 46 Only for a metal with 
relativistic dispersion, such as WSM and DSM, 
neutral Gaussian impurities can cause a R-linear 
LMC. In the presence of both ionic and short 
range scattering, the contributions to the LMR 
add up according to the the Mathiessen’s rule. 
Consequently, after entering the quantum limit 
(B > Bq), the LMR for a conventional metal with 
nonrelativistic dispersion first becomes negative 
and passes through a minimum at a nonuniversal 
field strength B = B*, before its eventual rise to 
become positive, as illustrated generically in Fig.Q] 

3. We explicitly discuss the Landau level (LL) struc¬ 
tures of various physical systems such as a three 
dimensional electron gas, a quasi two dimensional 
metal, a WSM, and a DSM with an even as well 
as an odd number of flavors. For each of these 
systems we show the validity of our general results 
regarding the LMC in the quantum limit. 

4. For layered materials with a corrugated Fermi 
surface, we provide a full quantum mechanical 
treatment for establishing the existence of the 
Yamaji angle in a tilted magnetic held with 
respect to the interlayer direction. We explicitly 
show that for a very small deviation from the 
Yamaji angle, a single LL remains partially filled, 
and the system effectively realizes a situation 
similar to the one in the quantum limit. This 
happens only for the layered materials with a 



B/flo 

FIG. 1. An illustration of the LMR due to the combined ef¬ 
fects of neutral and ionic impurities [Eq. 1341) 1. The quantum 
limit is achieved at a magnetic field strength B = Bq. The 
initial decrease of LMR for B > Bo is governed by ionic scat¬ 
tering and attains its minimum value at a nonuniversal field 
strength B * [Eq. 1351) 1. For B > B*, the sharp increase of 
LMR is controlled by the neutral impurities. The choice of 
B* = 2Bo has been made only for an illustrative purpose. 
The actual value of the crossover field B* depends on many 
nonuniversal sample-dependent details such as the compara¬ 
tive amounts of ionic long-range and neutral short-range dis¬ 
order, the carrier density, etc., and B * could, in principle, be 
large or small, thus suppressing the positive or negative LMR 
regimes in experiments. 

large carrier density, such that the energy scale 
associated with the areal density is larger than the 
interlayer hopping strength. Note that only for 
such a high carrier density, the underlying Fermi 
surface can be a corrugated cylinder, as found in 
PdCoC >2 and PtCo02- Other layered materials 
such as graphite^ and a-(BEDT-TTF) 2 l 3 can also 
exhibit a negative LMR in certain range of the 
field strength. However, these are compensated 
semimetals with very low carrier density and 
a set of closed electron and hole pockets. Our 
theory provides a simple explanation why the 
Yamaji angle and its associated phenomena will 
be absent in such materials with low carrier density. 

5. We discuss the experiments on Bii_ x Sb^, and 
ZrTeg^ in light of our theory. We show that the 
experimental observation of positive LMC oc B 2 
in these semiconducting materials are consistent 
with our theory in the quantum limit, with ionic 
scattering being the dominant relaxation mecha¬ 
nism, as is generically expected in semiconductors. 
We also argue why the observation of the positive 
LMC in ZrTeg^&, can not by itself, be construed as 
evidence for the existence of the chiral magnetic 
current. 
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The manuscript is organized as follows. In Sec. [XT] we 
demonstrate the relation between the axial anomaly as 
described by Eq. © and the electrical conductivity of a 
one dimensional system. We also show why the notion 
of the axial anomaly can be applied even for a bounded 
dispersion relation. In Sec. EU we explain how the one 
dimensional anomaly of a dispersive LL is responsible 
for giving rise to the axial anomaly (see Eq. ©) for a 
generic three dimensional system in parallel electric and 
magnetic fields. We demonstrate how the ABJ anomaly 
gives rise to the LMC and the LMR in contrast to the pre¬ 
dictions of the classical theory of magnetotransport. We 
again ensure that our conclusions are valid for bounded 
dispersion relations. In Sec. |TV] we calculate the LMC 
in the quantum limit, using a Boltzmann equation in the 
presence of short range and long range scattering due to 
neutral and ionic impurities respectively. In the Sec. |Vj 
we consider the LMC in the quantum limit for a three 
dimensional electron gas, as found in many semiconduc¬ 
tors. The discussion of the LMC for quasi-two dimen¬ 
sional or layered metals are presented in Sec. IVH We 
consider many aspects of the magnetotransport experi¬ 
ments on PdCo0 2 in Subsec. IVI Al The agreement be¬ 
tween our theory and the experimentally observed LMR 
in a-(BEDT-TTF) 2 l 3 is discussed in Subsec. IVI Bl The 
LMC of a WSM and a DSM with an even number of fla¬ 
vors are considered in Secs. IVIll and IVIIIl respectivelv. A 
brief discussion of the magnetotransport experiments on 
Cd 3 As 2 can be found in Sec. IVIIII Sec. [IX] is devoted to 
the discussion of a DSM with an odd number of flavors. 
In particular, we discuss recent magnetotransport mea¬ 
surements on Bii_ x Sb x and ZrTes in Subsec. IIX Al and 
Subsec. HX Bl respectively. In Subsec. UXBI we show why 
the observation of a positive LMC varying as B 2 does 
not necessarily imply the existence of a chiral magnetic 
current. Our main findings are summarized in Sec. [X] 
An explicit derivation of the transport lifetime using the 
Boltzmann equation is relegated to the Appendix [A] 


II. AXIAL ANOMALY OF ONE DIMENSIONAL 
METALS 

Quite amazingly, the axial anomaly described by 
Eq. © can naturally arise in any one dimensional elec¬ 
tron gas due to the existence of Fermi points and linearly 
dispersing low energy excitations in their vicinity. This 
one dimensional form of the axial anomaly will play a 
crucial role in the subsequent development of the theory 
of LMR in three dimensions in the presence of an external 
magnetic field. For this reason, we will first consider the 
implications of the axial anomaly for the electrical con¬ 
ductivity of a one dimensional system of spinless, charged 
fermions with dispersion e(fc). 

The low energy physics is governed by the linearly dis¬ 
persing right and left movers in the vicinity of the two 
Fermi points located at k = Ekp = ± 7 rn, where n is the 
linear density of fermions. The effective action for the 


linearized theory is given by 


S = 


dxdt [ihdt + ihvFTsd x \ ip- 


( 3 ) 


Here T 3 = diag(l, —1) is a Pauli matrix, and ip T = ( R , L) 
is a two component Grassmann spinor, with R and L re¬ 
spectively describing the right and the left movers, which 
are one component Weyl fermions. The Fermi velocity 
is defined as Hvf = dke(k)\k=k F - The linearized the¬ 
ory is invariant under two global gauge transformations 
ip -A e la ip and if> -A e iaaxT3 ip, which implies separate 
number conservation laws for the right and left movers. 
The total or electrical charge and current density opera¬ 
tors are defined as 

p{x,t) = —eip'ip, j(x,t) = —evF^Ts'ip. 

( 4 ) 


In contrast, the axial charge and current density opera¬ 
tors are defined as 


Pax (*G — 


JpM) 
evp ’ 


„• _ v F p(x,t) 

Jax — 


( 5 ) 


These relations between the total and the axial current 
operators are special properties of a two component Dirac 
fermion in one dimension. In the presence of an electric 
field E, the electrical charge and current still satisfy the 
continuity equation 


dtp + d x j = 0 . ( 6 ) 


In contrast, the axial current conservation law is violated 
according to Eq. ©, with Nf = 1. Notice that the spin 
degree of freedom can be taken into account by setting 
Nf = 2. In realistic solid state systems, Nf could be 
a larger integer because of valley degeneracy and band 
structure effects. We also note that the axial anomaly 
equation serves as the basis for the bosonization proce¬ 
dure. 

The physical consequence of the axial anomaly can 
be better understood by considering a constant external 
electric field. Such a field does not induce an electrical 
charge density p , which naturally sets j ax = 0. After in¬ 
tegrating both sides of Eq. © with respect to time, we 
find a linear growth of the axial charge and the electrical 
current with time, described by 


Pax — 


eVF 


eEt 
nh ’ 


( 7 ) 


for Nf = 1. The linear time dependence of j corresponds 
to a uniform acceleration of the center of mass. How¬ 
ever, in the presence of impurity scattering or a periodic 
potential due to the underlying lattice, this uniform ac¬ 
celeration can not be sustained for an indefinitely long 
time. For a generic solid state system, the transport life¬ 
time t tr due to the impurity scattering is much shorter 
than the Bloch oscillation period T = H/{eEa) due to 
the underlying periodic potential, where a is the lattice 
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spacing. Therefore, the anomaly equation will be perti¬ 
nent for a time scale t < r tr -C T. Since an axial charge 
build up is occurring between the two Fermi points, it is 
natural to anticipate the relaxation rate to be related to 
the backscattering rate. Based on this intuition, we can 
use the following phenomenological equation to describe 
the relaxation of the axial charge density 


&tPax 


p ax 
7 itl Tb 


( 8 ) 


In the steady state, this gives rise to a remarkable relation 
between the axial charge and the electrical conductivity 


_ eVFPax _ e 2 v F T b 
E 7 tH 


(9) 


where t^ 1 is the backscattering rate. We will be using 
this relation for the subsequent development of the axial 
anomaly for three dimensional metals and their magneto¬ 
conductivity. Therefore, we will first confirm the validity 
of this phenomenological description, by explicitly using 
the solution of the Boltzmann equation for the micro¬ 
scopic model, without linearizing the dispersion relation. 
This will also provide us with a concrete formula for the 
relaxation rate. 

Within the relaxation time approximation, the steady 
state current for the microscopic model is given by 


j = a E = 


e 2 T tr E f dk 


27r 


v 2 (k) 5(e F — e(A;)) 


e V F T tr 
irh 


E, 


( 10 ) 


where v(k) = dkt{k)/H is the group velocity. Notice that 
the same expression for the conductivity is predicted by 
the phenomenological description of the axial charge re¬ 
laxation if we identify Tt r = T b . Further confirmation of 
this comes from the solution of the linearized Boltzmann 
equation 


1 

Ttr 




v{k + q) 
v{k) 


( 11 ) 


The collision processes are captured by 

Ojr 

W(k,k + q) = —m\Uid{q)\ 2 S(e(k) - e[k + q)), (12) 
n 

where rii is the linear density of the impurities and U\d{q) 
is the Fourier transform of the one dimensional impurity 
potential. The linearized Boltzmann operator is anni¬ 
hilated by the forward scattering processes at q = 0, 
and the transport lifetime is entirely determined by the 
backscattering at q = —2k. This now provides a proper 
justification for identifying the transport lifetime with 
the backscattering lifetime, and we obtain the following 
expression 


h 2 v F 

Ttr =Tb= 2rii\Ui d (2k F )\ 2 " 


(13) 


Therefore, we can use the phenomenological descrip¬ 
tion of Eq. m, together with Eq. m , to describe the 
anomaly induced electrical transport in one dimension. 
Now we are in a position to address the relevance of ABJ 
anomaly for a three dimensional system in the presence 
of parallel electric and magnetic fields. 



FIG. 2. (Color online) The dispersing LLs of a three dimen¬ 
sional electron gas, where oj c = eB/m is the cyclotron fre¬ 
quency and Ib = y/h/eB is the magnetic length. The LLs 
with N = 0,1, 2 and 3 are respectively represented with the 
colors red, black, blue and green. The dashed brown lines cor¬ 
respond to two possible locations of the Fermi energy. When 
the Fermi energy lies entirely within the N = 0 LLL, we mark 
it by QL to illustrate that we are in the quantum limit. In 
the quantum limit only the partially occupied LLL produces 
a set of Fermi points. In contrast, away from the quantum 
limit, several partially filled LLs can produce Fermi points. 
Each partially occupied LL gives rise to the axial anomaly 
and participates in the longitudinal magnetotransport. 


III. AXIAL ANOMALY OF THREE 
DIMENSIONAL METALS 

We will first consider a generic three dimensional metal 
in the presence of a strong magnetic field B = (0,0, B). 
For simplicity we will assume the existence of a closed 
Fermi surface contour in the plane perpendicular to B. In 
this case, the cyclotron motion is quantized into discrete 
LLs, which can only disperse along the direction of the 
magnetic field. In the Landau gauge, A = (0,Ba;,0), we 
can describe the dispersion relation of the LLs according 
to 


e N {k y ,k z ) = e N (k z ), (14) 

which at least has the degeneracy eB/h , with a mag¬ 
netic length l F = \Jh/eB. There may be additional de¬ 
generacy of the LLs, arising from the presence of mul¬ 
tiple atoms in a unit cell, which we will denote by g n 
for each level. A partially filled LL at least crosses the 
Fermi energy at a set of Fermi points ±.k F ^(B) : and 
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acts as a highly degenerate quasi one dimensional con¬ 
duction channel. For example, we show a few LLs of a 
three dimensional electron gas in Fig. [2] When an exter¬ 
nal electric field is applied parallel to the magnetic field, 
there is no Lorentz force in this direction, and classically 
one does not expect any LMC or LMR. Only the par¬ 
tially filled LLs possessing Fermi points can participate 
in the longitudinal magnetotransport, and following the 
discussion in the previous section, we anticipate an im¬ 
portant connection between the one dimensional form of 
the axial anomaly and the LMC. 

Now consider the quantum limit, achieved for a very 
strong magnetic field strength B > B 0 , such that all the 
carriers are confined to the LLL with N = 0. The field 
scale Bq is determined by the density of carriers and the 
parameters of the underlying microscopic model. In the 
absence of impurity scattering, the axial anomaly formula 
in the quantum limit is obtained by multiplying the right 
hand side of Eq. m by eB/h and identifying Nf as the 
orbital or spin degeneracy go in the LLL, which leads to 


e 2 EB 

OtPaxfi ~ 50 27r2ft2 - 


(15) 


Notice that we have obtained the ABJ anomaly formula 
of Eq. ([2]) without referring to any underlying relativistic 
band structure. Therefore, the emergence of the axial 
anomaly is a universal feature of all three dimensional 
metals placed in parallel electric and magnetic fields. In 
the presence of impurity scattering, the axial charge and 
the magnetoconductivity in the quantum limit can be 
determined by following Eq. ©. They are given by 


Pax, 0 — 90~ 


B) = go 


e 2 EBTb t o(B) 

e 2 VF,o(B)T b ' 0 (B) 
27 T 2 hl 2 B 


(16) 

(17) 


where r b q 1 is the backscattering rate between the Fermi 
points and v F ,o (B) is the Fermi velocity in the LLL. 
The backscattering rate, which depends on the nature 
of disorder, has to be determined from the Boltzmann 
equation for longitudinal magnetotransport in the LLL. 
The evaluation of t^o is shown in the following section. 
The formula for LMC in Eq. (fT71) is one of the main re¬ 
sults of this manuscript and is applicable to any weakly 
correlated generic metal in the quantum limit, indepen¬ 
dent of the explicit form of the band dispersion. This is 
the generalization of Nielsen-Ninomiya’s formula^ for ax¬ 
ial anomaly induced LMC of Weyl fermions to a generic 
three dimensional meta. As long as v Fy o{B)T b (B)/l 2 B is 
not a constant, the axial anomaly can give rise to LMC 
and LMR, in contrast to the anticipations of the classi¬ 
cal theory of magnetotransport. The issue of positive or 
negative LMR is determined by the explicit dependence 
of the relaxation rate on B , depending explicitly on the 
nature of the underlying disorder, which we will consider 
in the next section. 

It has been known for a long time that the dimen¬ 
sional reduction caused by the strong magnetic field in 


the quantum limit, can also lead to a charge density wave 
(CDW) instability of the electronic system at very low 
temperatures^ - —. Such instabilities generically tend to 
gap out the Fermi points, and will lead to an insulat¬ 
ing behavior of the LMR for a clean, interacting system. 
However, it is important to note that an experimental 
detection of such a broken symmetry phase is quite chal¬ 
lenging and signatures of certain instabilities have only 
been reported in graphite and bismuth. Our theory is 
based on the existence of the Fermi points in the absence 
of any such instability. This assumption can be justi¬ 
fied in the following way. The disorder potential couples 
to CDW order as a random field, and it generally plays 
a detrimental role for the existence of a CDW as a true 
long range order. If we consider a continuum model of the 
electron gas, a charge density wave breaks the continu¬ 
ous translational symmetry, and the Imry Ma argument 
forbids the existence of a continuous symmetry break¬ 
ing below four dimensions in the presence of a random 
field. In contrast, for a lattice model of electrons, a CDW 
only breaks a discrete translational symmetry, and the 
Imry Ma argument does not prohibit its existence as a 
long range order in three dimensions. However, disorder 
can cause a severe suppression of the transition temper¬ 
ature^, while making it effectively impossible to detect 
such ordering in most materials. In any case, all our re¬ 
sults would remain valid at temperatures higher than the 
transition temperature of such a CDW instability (in the 
unlikely scenario that an instability does occur), and we 
ignore the possibility of a CDW instability in the current 
work. 

For weaker magnetic field strengths below the quantum 
limit B < So, there are multiple partially occupied LLs. 
Each of these levels will have their own ABJ anomaly 
defined by the equation 


„ _ e 2 EB 

OtPax,N ~ gN 2?r 2^2 ’ 


(18) 


where gN is the orbital and/or spin degeneracy in 
the ./Vth LL. Thus each LL contributes to the LMC 
through its transport lifetime T tr ,N{B) and Fermi veloc¬ 
ity v f ^n{B). This leads to a very general formula for the 
LMC ’ 


T (B) =J> 

N 


e 2 V F , N (B)T tr ,N(B) 


2n 2 hl 2 B 


(19) 


where the prime denotes a sum carried over all the par¬ 
tially filled LLs. The transport lifetimes of various LLs 
follow a set of coupled linear algebraic equations, which 
account for both intra and inter LL scattering processes. 
This will be derived in the Appendix [A] In the Boltz¬ 
mann theory, the broadening of the LLs is not taken 
into account. Therefore, Eq. (fl9l) will lead to undamped 
Shubnikov de Hass oscillations of the LMR. The inclu¬ 
sion of the broadening at a phenomenological level will 
have two important consequences: (i) a Dingle factor 
for each level, which will suppress the amplitude of the 
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oscillations, and (ii) a quantitative modification of the 
transport lifetimes due to the forward scattering process. 
However, the backscattering processes with large momen¬ 
tum transfer still make the dominant contribution to the 
transport lifetime. However, these quantitative details 
do not alter the fact that the axial anomaly gives rise to 
the very existence of a finite LMR. 

We can further demonstrate that Eq. (USD is not an 
artefact of a linearized theory, and is valid for an arbi¬ 
trary LL dispersion. In the presence of parallel electric 
and magnetic fields, the current carried by the quasipar¬ 
ticles residing in the partially filled LLs are given by 

■’ z = ~2^? / ~^r^29N VN ( k z)f N {k z ). (20) 

In the relaxation rate approximation, the deviation from 
equilibrium is described by 


Of o 

fN(k z ) = f N {kz) - eEv N (k z )T tr , N {k z )-^-\ e=eNikz) . 


After substituting this in Eq. (EU1) . we obtain 


( 21 ) 


Jz = 


e 2 E 
2n 2 hl 2 B 


9NVF,NT t tr,N (B ), 
N 


( 22 ) 


which indeed leads to the most general expression for the 
LMC announced in Eq. ©■ 

It is possible to reach the quantum limit without using 
very large magnetic fields in certain narrow-gap and gap¬ 
less semiconductors with very low carrier densities. It is 
also possible to attain the condition of a single partially 
filled LL in a quasi-two dimensional metal with high car¬ 
rier density. This condition can be realized either by 
making the cyclotron energy bigger than the interlayer 
hopping strength or adjusting the angle between the ap¬ 
plied magnetic and electric fields in the proximity of the 
“Yamaji angles”. We discuss these aspects in detail in 
the Sec. I VI I In the next section we focus on determin¬ 
ing the transport lifetime and the LMC in the quantum 
limit. 


IV. LONGITUDINAL 
MAGNETOCONDUCTIVITY IN THE 
QUANTUM LIMIT 


Computing the transport lifetime from the Boltzmann 
equation, within the relaxation time approximation for 
a generic dispersion £iv = o(fc z ) in the quantum limit, we 
arrive at 


_ k 2 VF,o(B) 

bA j 2n l l 2 B \U eff {2k Ffi m\ 2: 1 j 

where n* is the impurity density per unit volume and 
U e ff (q z ) is an effective one dimensional scattering po¬ 
tential. This scattering potential has to be obtained 


by integrating the matrix elements of the three dimen¬ 
sional disorder potential over the in-plane components of 
the momentum, and contains information regarding the 
wavefunction of the LLL. Remarkably, the result of the 
Boltzmann calculation in the presence of the magnetic 
field agrees with the one dimensional result of Eq. dED, 
where the one dimensional potential U\d is replaced by 
U e f f. The expression for the effective potential is given 

by 


I U ef f{q z )\ 2 = J j^^\U(q ± ,q z )\ 2 exp , 

(2d) 

where U(q±,q z ) is the Fourier transform of the actual 
three dimensional real-space disorder potential. We now 
consider two types of disorder potentials: (i) short range 
Gaussian impurity scattering specified by U s (q± : q z ) = 
Uq exp[— a 2 (qj_ + q 2 )/ 2], where a is the range of the im¬ 
purity potential, and (ii) long range ionic impurity scat¬ 
tering characterized by the screened Coulomb potential 
U c (q±,q z ) = (47re 2 / K)[q\ + q 2 + < 7 y F ] -1 , where k is the 
background lattice dielectric constant of the system and 
1/qTF is the Thomas-Fermi screening length. We note 
that in the limit a —> 0 the Gaussian potential reduces to 
an onsite Dirac delta function potential UoS(r) represen¬ 
tative of the zero-range disorder associated with neutral 
point defects. 

For the Gaussian potential, we find the following effec¬ 
tive one dimensional potential and the transport lifetime 


U!ff(Qz) = 


Uo 


\J2'k1 2 b {1 2 b 


2a?) 


exp - 


2 2 
& 


ttE 


rf,o( B ) = - 7 yj2V Ffi (B)(l B +2a )exp[4 k Ffi (B)a 

U n 


(25) 


(26) 


where nf is the density of neutral Gaussian impurities. 
Therefore, the contribution of the Gaussian short-range 
impurity scattering to the LMC in the quantum limit is 
given by 

p 2 h f 2 n 2 \ 

= ^u! v Ffi{B) [1 + jr) exp[4 k 2 Ffi {B)a\ 

(27) 


Depending on the underlying dispersion relation, vf,o can 
be a periodic function of kp,o(B). This in turn can cause 
a nonmonotonic dependence of cr{B) on the magnetic 
field strength. A detailed discussion of the nonmono¬ 
tonic variation of the LMC due to zero range impurity 
scattering or point defects (a —> 0) will be provided in 
the following sections, when we consider different phys¬ 
ical systems. However, for a system with nonrelativis- 
tic dispersion and very large magnetic field strength, the 
asymptotic behavior of <J S {B) is given by 

a s (B)~l 2 B (l 2 B +2a 2 ). 


Therefore, the short range scattering due to neutral im¬ 
purities always lead to positive LMR for a conventional 
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FIG. 3. LMC for an electron gas in semiconductors [Eq. (1371) 1 in the quantum limit, (a) Gaussian potential scattering: LMC 
with a fixed carrier density of n = 6.75 x 10 26 mA 3 as a function of (o/fs) 4 , where a is the range of the scattering potential, 
and <r s (0) is the conductivity in the absence of a magnetic field due to short range scattering (see the Appendix). We find 
a s (B)/a s ( 0) decreases as a function of B, which gives rise to a positive LMR. (b) Ionic scattering: LMC with a fixed carrier 
density of n = 6.75 x 10 26 m -3 as a function of (oo /Ib) 4 , where a o is the Bohr radius in the vacuum and ct c ( 0) is the conductivity 
in the absence of a magnetic field (see the Appendix). We find a c (B)/a c (0 ) ~ B 2 , which gives rise to a negative LMR. 


metal in the very large magnetic field limit B —> oo. In 
contrast, the finite range Gaussian impurities can lead 
to a H-linear positive LMC for systems with underlying 
relativistic dispersion, which will be discussed in the fol¬ 
lowing sections. 

For long-range (i.e. Coulombic) ionic impurity scatter¬ 
ing, the effective one dimensional potential is given by 

U c eff {q z ) = ^^Jr[-lAql + q 2 TFi0 (B))ll/2} X 

exp [(Qz + QtfA B )) * 1 b/ a P 8 ) 

where r(— 1, x) is the upper incomplete Gamma function, 
and q FF 0 is the Thomas Fermi screening length in the 
quantum limit. By using this effective potential we find 
the following contribution due to the ionic impurities to 
the LMC 


Hti 2 v 2 Ffi (B) exp 

2 {4^F, q{B) + Qtf, o(-®)} 

8^ 3 e 2 n*l% r 


1,|{4 kl 0 (B) + q 2 TFfi (B)} 


(29) 

where n\ is the density of the charged impurities. The 
analysis of the asymptotic behavior of <j(B) can be facil¬ 
itated by using the following relation 

OO 

e x T{-l,x) = ^](-l) n r(n + 2)^- (n+2) . (30) 

n =0 . 


In the quantum limit, the enhanced degeneracy of the 
LLL leads to very strong screening of the ionic impurities, 
which is captured through the following field dependent 
Thomas Fermi momentum 


qrF,o(B) 


e 

^/2tt 2 Hv F}0 (B)1 2 b 


(31) 


Since q FF ,o an d kp ,o are respectively increasing and de¬ 
creasing functions of B , the asymptotic expansion is dom¬ 
inated by qrF, o- This leads to a strong positive LMC in¬ 
creasing with field strength as ~ B 2 . This is the generic 
high-field, extreme quantum limit behavior of charged 
disorder induced LMC in three-dimensional metals. 

However, this asymptotic behavior can be understood 
on a more physical ground. Since the screening length 
in the quantum limit decreases with the increasing mag¬ 
netic field strength, the effective one dimensional scatter¬ 
ing potential also decreases as U ^ ~ ^tfo^b 1 - Conse¬ 
quently, the relaxation time t£ 0 oc v F ,o ( B eff)~ 2 increases 
with B. As the conductivity in the asymptotic limit 
B —> oo behaves as <J C (B) ~ ^fo9tfoj using Eq. (13T1) 
we arrive at one of our key results 

a c (B) ~ l~ B A ~ B\ 


which is valid for an arbitrary dispersion relation of the 
LLL. When ionic scattering is the only resistive relax¬ 
ation mechanism, there is a large negative LMR in the 
quantum limit, and the positive LMC u c (B) oc B 2 de¬ 
scribes a universal result. 


is given by Matthiessen’s rule 

1 _ 1 , 1 

n,o(B) ~ t£ 0 (B) + T § fi (B)' 


In the presence of both types (i.e. long-range ionic and 
short-range neutral) of impurities, the net relaxation rate 


( 32 ) 
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which leads to the net LMR 


p(B)=p c (B)+p s (B). (33) 


Here p c (B) = 1 /a c (B) and p s {B) = 1 /cr s (B). When 
a c (B) < a s (B ), we obtain a negative LMR and a posi¬ 
tive LMC cr(B) oc B 2 . For magnetic field strengths above 
a nonuniversal threshold B *, the contribution from short 
range impurity dominates, which leads to a positive LMR 
for a conventional metal. Therefore, short range impuri¬ 
ties always become important for an asymptotically large 
field strength. The value of B* depends on the ratio of 
the impurity densities (n£/nf), and various material pa¬ 
rameters such as effective mass and density of carriers. 
For concreteness, if we consider only zero-range point im¬ 
purities as the source of short range scattering, we find 


P{B) 


p{B*) (B *\ 2 

2 



h 2 y/ivn f n- \ 

ey/U 0 Km z \n*J 


(34) 

(35) 


where n is the density of the carriers and m~ is the effec¬ 
tive band mass along the direction of the applied field. 
Notice that we are defining B* through the condition 
p c (B*) = p s (B*) = p(B*)/2. The generic behavior of 
the LMR described by Eq. (I3H) is illustrated in Fig. [Q 
In the following sections we establish the applicability 
of our theory for various experimentally relevant mate¬ 
rials. We consider different physical systems according 
to the varying degree of complexity of the underlying 
band structure. We begin with the simplest situation of 
a three dimensional electron gas with parabolic disper¬ 
sion, followed by the consideration of a simple quasi two 
dimensional metal. In the absence of the magnetic field, 
both systems have single sheet of Fermi surface. Later we 
consider the examples of WSM and DSM, which can pos¬ 
sess multiple Fermi surfaces due to the additional orbital 
degrees of freedom and spin-orbit coupling. 


V. ELECTRON GAS IN SEMICONDUCTORS 

Usually semiconductors have low carrier density, and 
for this reason it is easier to achieve the quantum limit 
even with a moderately strong magnetic field. Within the 
k ■ p theory, the carriers in many semiconductors can be 
modeled by a conventional three dimensional electron gas 
with a parabolic dispersion. The pertinent Hamiltonian 
in the presence of the magnetic field is (ignoring spin) 

Hi = -^-(p-eA) 2 , (36) 

2m 

where m is the effective band mass. In this case, the LLs 
have the following dispersion relations 

/ 1\ h 2 h 2 

e N (k y ,k z )=(N+-\hu c +y^, (37) 


where the cyclotron frequency co c = eB/m. First four 
dispersing LLs and their possible intersections with the 
Fermi level are illustrated in Fig. [2] The quantum limit is 
achieved when the field strength B exceeds the threshold 
value 

B a = - (27tV) 1/3 . (38) 

The Fermi wavenumber for this problem decreases with 
B according to 

171V fo „ 2 i2 /ore 

k F ,o = — z — = 2tt nl B , (39) 

n 

and the Thomas Fermi screening length is determined by 


-l _ 

Qtf, o — 



(40) 


For short range impurity scattering, we find 


2e 2 TT 3 a 2 n 2 h 3 2 
a{B) = n.t/fr.i 


1 + 


/ 2 

l b 

2 a 2 


exp(167r 4 rrGrZ]j), 


(41) 

which is a monotonically decreasing function of B (see 
Fig. 3(a)). For a finite range of the impurity potential, 
the LMC deceases as ~ B~ 3 . While for zero-range point 
impurities (a —>■ 0), the LMC decreases as ~ B~ 2 . In 
contrast, the axial anomaly in the presence of long-range 
ionic impurity scattering l eads to a large positive LMC 
growing as ~ B 2 (see Fig. |3(b)| ). After substituting the 
expressions for k_ F ,o, v F p and qr F , o hr Eq. (12911 . we obtain 
an explicit form of the magnetoconductivity due to the 
ionic impurities, which precisely agrees with the formula 
derived by Adam and Argyres^. After combining the 
contributions from both types of impurity scattering, we 
obtain the net LMR given by Eq. where 


B* = Bn 


1.92 h 




o Km n 


1/3 


1/4 


(42) 


Therefore, a low density of carriers simultaneously de¬ 
creases Bq and increases the crossover field strength B*, 
which is suitable for observing the negative LMR over 
a significant range of the applied field strength. Due to 
the presence of donors and acceptors in semiconductors, 
the scattering by ionic impurities is the most important 
relaxation mechanism at low temperatures. Therefore, 
semiconducting systems will provide many suitable can¬ 
didate materials for observing the axial anomaly induced 
negative LMR in the quantum limit. In addition, semi¬ 
conductors, because of their low carrier density, can be 
easily driven into the quantum limit of the LLL occu¬ 
pancy, necessary for the generic axial anomaly being dis¬ 
cussed in this work. 

An explicit consideration of the Zeeman spin splitting 
does not change our conclusions. In the absence of spin 
flip scattering, there is no impurity induced matrix ele¬ 
ments between two Zeeman split LLs in the LMC cal¬ 
culation. Therefore, our calculations go through for each 
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FIG. 4. (Color online)The dispersive LLs for a quasi-two dimensional metal [Eq. (1461) 1. when the magnetic field is applied along 
the c axis, (a) 4 t z > hw c and (b) 4t z < hw c . For a high carrier density corresponding to Eq. (1441) . the underlying Fermi surface 
is a corrugated cylinder and hw c exceeds 4 t z before reaching the quantum limit. In this case, each LL is well separated as 
shown in panel (b) which gives rise to two Fermi points, and the system behaves effectively like the one in the quantum limit. 


spin split LL without any essential modification. We only 
need to ensure that the total density of carriers is dis¬ 
tributed between two spin split levels. In the very strong 
field limit, eventually one arrives at a fully spin polarized 

I 


VI. QUASI-TWO DIMENSIONAL METALS 


In this subsection we consider the axial anomaly in¬ 
duced LMC of quasi-two dimensional metals. By quasi- 
two dimensional metals we are referring to layered ma¬ 
terials where the inter-plane hopping matrix element is 
much smaller than its in-plane counterparts. Many in¬ 
teresting materials such PdCo 02 , PtCo 02 , graphite and 
a-(BEDT) 2 l 3 belong to this category. We will focus our 
discussion in Subsec. IVI Al on PdCoC> 2 , which has a single 
sheet of Fermi surface^ — —. In Subsec. IVIBI we con¬ 
sider the qualitative agreement between our theory and 
the longitudinal magnetotransport measurements on a- 
(BEDT) 2 l 3 - Towards the end we also provide a brief 
discussion of experiments on graphite. 


A. PdCoCU and PtCoCU 


For simplicity we consider the following Hamiltonian 

Hq2D = 7T^- 2tz cos(k*c) 3 (43) 

2m± 

where p± = (p x ,Py, 0), m± is the effective mass in the 
ab plane, and c is the lattice spacing in the interlayer 
direction. When interlayer coupling is sufficiently weak 
in comparison to the energy scale of the in plane motion 


limit, which is exactly described by our calculations for 
a spinless model. Therefore, our results properly capture 
the LMC in the asymptotic limit B —> oo. 


e., 


tz < 


nH 2 nc 

2m± 


(44) 


the underlying Fermi surface is a corrugated cylinder. 
This is the situation for PdCoC> 2 , PtCo02. In contrast 
for a stronger interlayer hopping such that 


tz > 


7r H 2 nc 
2m± ’ 


(45) 


the system possesses an ellipsoidal Fermi surface. In this 
case, one can expand the cos (k z c) and obtain an effective 
parabolic dispersion like the one for an electron gas. Even 
though graphite has a more complicated band structure 
due to the presence of multiple atoms in a unit cell, which 
gives rise to a compensated semimetal phase with both 
electron and hole pocket s 54 ’ 55 , some of the qualitative 
features can still be understood in terms of an electron 
gas with an ellipsoidal Fermi surface. In particular there 
are some interesting differences between PdCoCU and 
graphite in a tilted magnetic field, which we can easily 
understand on the basis of a simple Hamiltonian H Q 2 d- 
When the magnetic field is applied along the c axis, 
the dispersion relations for the LLs are described by 


CN ij^'y 5 



hu> c — 2t z cos (k z c), 


(46) 


where uj c = e /i /m ±. T he dispersive LLs for 4f 2 > Huj c , 
are shown in Fig. 4(a)| when the Fermi level can cross 
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FIG. 5. (Color online)The dispersive LLs of a quasi-two dimensional metal with a high carrier density in a magnetic field 
(solid) applied at an angle 8 with respect to the c axis, as a function of k' z c (defined in the text), for t z /(huj c ) = 2. We are 
considering a large Fermi energy such that the Fermi level crosses only the LLs with large index N. For (a) 8 = 0 there are a 
large number of Fermi points from multiple LLs. (b) For the same field strength and tilt angle 8 = 0.097r the LLs with high 
index N become completely flat and do not produce any Fermi point. Therefore the effective hopping strength along the field 
direction, for a large LL index vanishes, showing the existence of the Yamaji angle, independent of the magnitude of t z /(hui c ). 
(c) For a small deviation away from the Yamaji angle, at 8 = 0.087T there is only one partially filled LL producing two Fermi 
points, similar to the situation in the quantum limit. 


multiple LLs similar to the situation in an electron gas. 
In contrast, when 4 t z < Hu} Cl all the LLs are well sepa¬ 
rated, and the Fermi energy can intersect only one LL at 
a time, which is illustrated in Fig. 4(b) For a material 
like PdCo02, which satisfies Eq. (1441) . Tujjn has to exceed 
the scale 4f, before the magnetic field can confine all the 
carriers into the LLL (true quantum limit). This is a 
unique situation when one almost mimics the quantum 
limit for a smaller field strength, where the axial anomaly 
due to a single LL contributes to the LMC. In this case, 
the inter LL scattering effects do not appear in the Boltz¬ 
mann calculation (similar to the quantum limit), and we 
only need to use the effective potential 


\U e ff,N{qz)\ 2 = J \U(q±,q z )\ 2 exp 


(2tt) 2 Z! 

r , q 2 J% 

Ln 




(47) 


for the corresponding partially occupied LL. In this case 
for short range scattering it is possible to find a closed 
form expression for an arbitrary LL. For ionic scatter¬ 
ing we can also obtain the corresponding expressions in 
closed form, but they become progressively lengthier with 
the LL index N. Some of these details can be found in 
the Appendix. When 


2'Kncl 2 B < 1 , 

I 


In a material like PdCo02, to attain At z < Huj c , one 
requires B > 96T4£ Therefore, multiple partially occu¬ 
pied LLs contribute to the LMC for any experimentally 
accessible field strength. The calculation of LMC using 
Eq. (TT!?!) is quite involved and it is not easy to unambigu¬ 


only the LLL is partially occupied, and the LMC in the 
quantum limit due to the short range scattering is given 
by 


a s (B) 


4e 2 f 2 c 2 . 2 2 ,2 \ ( 2 a 2 \ 

x exp[167r 2 ?z 2 /^a 2 ], 


(48) 


where we have used kp ,o = 2TT 2 nl B , and vf ,o = 
2t z csm(kF,oc)/h. Due to the periodic dependence of the 
Fermi velocity on B , we have a nonmonotonic behavior 
of this contribution to the LMC, which initially increases 
with the magnetic field strength until 47 rncl B < 1. For 
a stronger magnetic field, the LMC due to short range 
scattering decreases. When, Airndg -C 1, the sin 2 term 
can be approximated as ~ l B , and produces the same 
asymptotic behavior as in the case of an electron gas. 
Therefore, in the presence of only short range scatter¬ 
ing, one may initially observe a negative LMR, followed 
by an upturn to a positive LMR, in the quantum limit. 
In contrast, the ionic scattering always leads to a strong 
positive LMC varying as B 2 . Therefore, in the presence 
of both types of impurity scattering, the axial anomaly 
induced LMR follows the generic behavior illustrated in 
Fig. Q] 


ously relate the observed negative LMR with the axial 
anomaly at a quantitative level. How can we then probe 
the large negative LMR due to the axial anomaly of a 
single LL? This can be achieved by placing the system in 
a tilted magnetic field. Due to the high carrier density in 
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each layer, for any experimentally accessible steady field 
strength the Fermi level crosses only the LLs with large 
index N. When the field is tilted away from the c axis, 
the effective interlayer hopping is diminished, and for cer¬ 
tain tilt angles Oy , also known as the Yamaji angles, the 
LLs with large N become completely flat or nondisper- 
siv o 44 ’ 45 , which can not produce any Fermi point. A 
slight deviation from 0y can restore a weak dispersion 
of these LLs along the applied field direction. Conse¬ 
quently, even for a moderate field strength we can have 
a single partially occupied LL, similar to the situation in 
Fig. |4(b)| The original calculation of Yamaji^i is a semi- 
classical one, and does not explicitly consider the Fermi 
points produced by the dispersive LL’s, which play the 
most important role in our theory. Subsequently, the 
dispersion relations of the LLs have been calculated per- 
turbatively in the limits t z /(hoj c cos9) << 1, where 6 
is the angle between the c axis and the magnetic field 
B. However, the experiments are performed in the op¬ 
posite limits t z /fujj c >1. In addition, when 9 becomes 
very small, the control parameter always becomes much 
bigger than unity. For this reason, we have numerically 
solved the relevant Schrodinger equation with an arbi¬ 
trary value of t z /hu) c , from which we have established 
the existence of 9y as a nonperturbative quantum me¬ 
chanical effect. 

We consider a tilted magnetic field B = 
(0, B sin 9, B cos 9) making an angle 9 with the c axis and 
choose the vector potential A = (0,xB cos9, —xBsin9). 
After the Peierl’s substitution k —> k eA/h, we arrive 
at the following Hamiltonian 

h 2 2 h 2 2 

H q 2D = - o x + - (ky - mu> c cos 9x) 

Zm± Zuij_ 

— 2 t z cos(c[fc 2 + muj c sin 9x\ ) . (49) 


The momentum component along the field direction k' z = 
sin 9k y + cos 9k z and k y by itself are good quantum num¬ 
bers. This becomes transparent, when we transform the 
x coordinate to x' = x — k v / ( m±u> c cos 0), which leads to 


h 2 1 

Hq 2 D = —Z - d 2 , + -my(hx c cos9x'Y 


2m± x 2 
— 2t z cos(c' [k' z + mj_Lu c cos 9 sin 9x'} ), (50) 


where c' = c/ cos 9. We numerically obtain the eigenval¬ 
ues and the eigenstates of this Hamiltonian in Eq. (l50l) 
after discretizing x'. 

We consider a model parameter t z /(hui c ) = 2, which 


lies well outside the perturbative regime. The Fig. 5(a) 


corresponds to the field along c axis, where many LLs 
cross the Fermi energy. The Fig. |5(b)] corresponds to the 
field tilted at the first Yamaji angle 9 = 0.09-7T. At this 
angle, we clearly see the existence of non-dispersive LLs 
with high index N. For a small LL index N, the levels 
are still dispersive. This is the reason, why in contrast 
to the dense delafossite materials, a dilute system like 
graphite can not show the existence of Yamaji angle in 
AMRO measurements. It is important to note that there 


exist multiple Yamaji angles and within a perturbative 
treatment these are given by the zeros of a Bessel func- 
For a small deviation from the first Yamaji 


tio: 


,44,45 


angle, we find a single, weakly dispersive LL intersect¬ 
ing the Fermi energy, which is shown for 9 = 0.087T in 
Fig. 5(c) This yields only two Fermi points similar to 
the situation in the quantum limit. 

When the field is tilted, the measured magnetoresis¬ 
tance along the c axis is a combination of the LMR along 
the field direction and the TMR in the plane perpendicu¬ 
lar to B. At the Yamaji angle, there are no dispersive LL 
with Fermi points, and we do not have any axial anomaly. 
For this reason, at the Yamaji angles, the magnetoresis¬ 
tance along the c axis is completely determined by the 
positive TMR in the plane perpendicular to the applied 
B. In contrast, for a small deviation from the Yamaji an¬ 
gle, the axial anomaly of a single dispersive LL will give 
rise to a LMR along the B. In experiments on PdCoC> 2 , 
a small deviation from the 9y reveals the emergence of a 
negative component of the magnetoresistance measured 
along the interlayer direction. Therefore, axial anomaly 
in this material indeed leads to a negative LMR. More 
striking aspect is that for a sufficiently large B , the nega¬ 
tive LMR actually subdues the very large positive TMR. 
At present we do not have any quantitative understand¬ 
ing of this effect. But, the observation certainly provides 
evidence for a very large magnitude of the anomaly in¬ 
duced negative LMR from a single partially filled LL, as 
anticipated in the quantum limit. 


B. q-(BEDT-TTF) 2 I 3 

a-(BEDT-TTF) 2 l 3 is a layered material, which under 
high pressure realizes a semimetallic stated. The in¬ 
tralayer dispersion is captured by two tilted Dirac cones 
and the system is half-fille d 57 ^ 58 . Without the tilt, the 
in-plane dispersion is analogous to the quasiparticle spec¬ 
trum of graphene. Here the band degeneracy is acciden¬ 
tal, and do not occur at any high symmetry points. A 
negative LMR and a positive TMR have been concomi¬ 
tantly observed in this material for B > 0.7 T (Refs. [48l 
and 159). For understanding the salient features of the 
longitudinal magnetotransport measurements, we can ig¬ 
nore the tilt of the Dirac cone in the ab plane as in Ref.|59j. 
But, instead of a perturbative expansion of conductivity 
in interlayer hopping strength t z , we retain the full dis¬ 
persion along the c-axis. 

We consider the following effective Hamiltonian 

H or = hvk± ■ t ® r /3 — 2t z cos (k z c), (51) 

where the Pauli matrices Tj and rjj respectively act on the 
sublattice and the valley indices. The dispersion relations 
for the conduction and the valence bands are given by 

£k = ±hvk± — 2t z cos (k z c). (52) 

From the dispersion relation we can show that the density 
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of states behaves as 

d{€) “ W' (53> 

Therefore, for a system at half-filling the thermally ex¬ 
cited carrier density (both electron and holes) will be 
nr ~ T , in agreement with experiments. Due to the 
very low carrier density the conduction and the valence 
bands give rise to electron and hole pockets respectively 
centered around k z = 0 and k z = 7 r. At half-filling 
(ep = 0), the volume of the electron and hole pockets 
(elongated along the c axis) are equal, and they both 
end at k z = ±7 t/2 . Thus, the system behaves as a com¬ 
pensated semimetal. 

In the presence of a quantizing magnetic field, the LL 
dispersion around each valley is given by 

£iv#o,fc 2 , s = ±V2N ^ - 2 t z cos(fc 2 c) - s A Z: (54) 

lb 


the one in Fig. [H with the exception that the LMR finally 
goes to infinity for A z > 2 t z . 

Before leaving the topic of quasi two dimensional met¬ 
als, we will also briefly mention the experimental situa¬ 
tion in graphite. Due to the very low carrier density, it 
is easy to attain the quantum limit with a moderately 
strong magnetic field B ~ 7.5 T along the c axis— At 
a stronger magnetic field graphite suffers an instability 
towards a charge density wave formation with a wavevec- 
tor 2/cf,o(-B, n) around B ~ 25T (Refs. l47l and l49l). For 
9 T < B < B C (T), where B C (T) is the threshold for the 
first instability of the ground state, the magnetotrans¬ 
port data shows concomitant negative LMR and positive 
TMR (see Fig. 1 of Ref. d?]). This range becomes larger 
with increasing temperature. A more sophisticated treat¬ 
ment that accounts for electronic interaction is required 
for addressing the magnetotransport when B > B C (T ). 

VII. WEYL SEMIMETAL 


where s = ±1 denote the spin projections and A z = 
gfMsB/2. In contrast the dispersion for the zeroth LL is 
captured by 

ejv=o,fc*,s = -2 t z cos (k z c) - s A z . (55) 


When 


V2— > (4t + 2A Z ), 

IB 


(56) 


only the LLL can be partially filled. At half-filling, the 
spin down and the spin up LLL become completely empty 
and filled respectively if the Zeeman splitting exceeds 
the band width 2 t z , and the LMC vanishes (system be¬ 
haves as a band insulator). Therefore a finite LMC in the 
quantum limit can only occur for 2 t z > A z . The Fermi 
wavenumbers in the LLL are given by 


kp a =-b s — arcsin 

2c c 



(57) 


The LMC due to the short range scattering is now given 

by 


\B) = 


SeW 

hniU* 

(g4^F,-t- a 


1 + 


2a 2 

W 

l B J 


+ e 


iki 


1-^1 

4 £ 2 . 

)0(2tz - A z ), (58) 


where 0(:c) is the Heaviside step function. Therefore, 
zero range point impurities can only cause a positive 
LMR. In contrast, for weak enough magnetic fields the 
Gaussian impurities can lead to a I? linear positive LMC, 
and in this regime the negative LMR will vary as p ~ 
1 /(B* + B), which has been observed in Ref. |48|. In the 
presence of ionic impurities, the LMR can initially de¬ 
crease as 1/B 2 until the suppression due to the Zeeman 
splitting becomes strong enough to cause a positive LMR. 
Therefore, the overall profile for the LMR is similar to 


WSM is a gapless state where the nondegenerate con¬ 
duction and the valence bands touch linearly at isolated 
points in the Brillouin zone. For this reason, the touching 
points act as the monopole and the antimonopole of the 
U(l) Berry flux, and lead to many intriguing transport 
and optical properties^ - — . One must break time rever¬ 
sal or inversion symmetry for realizing a WSM. There are 
many theoretical proposals for realizing a WS M 13 ' ' 63 - —. 
But, a concrete experimental example of a noncentrosym- 
metric Weyl semimetal phase has only been found very 
recently in TaAsii—, and the band structure calcula¬ 
tions also predict the existence of such a phase in TaP, 
NbAs and Nb P 14 ’ 15 . There is also some preliminary ev¬ 
idence for a WSM in /?-Ag 2 Se— Given that a massless 
Dirac fermion is made up of two Weyl fermions of oppo¬ 
site chirality, which are located at the same point in the 
Brillouin zone, and a Zeeman coupling explicitly lifting 
the spin degeneracy can convert a DSM into a WSM, it is 
important to first understand the magnetotransport for 
the WSM. 

For this reason, we consider the following toy model 
on a tetragonal lattice 

Hwsm = 2£_L,i[sin(/c x a)<7i + sin(fej,a)o - 2 ] - [2£_l, 2 {2 - 

cos(fca,a) — cos(fc y a)} + 2 t z cos (k z c) — A]ct 3 , (59) 

to describe a WSM, where a and c are respectively the 
lattice spacings in the ab plane and along the c axis. For 
\A/t z \ < 1, conduction and valence bands touch at 

k= ^0,0, ±-arccos(A/2f 2 )^ , (60) 

which are the respective locations of the left and the right 
handed Weyl fermions. The tight binding parameter f j _,2 
prohibits the existence of additional gapless points for 
t / 0 and k y ^ 0. For studying the LLs of the Weyl 
fermions, we linearize F 2 in k x and k yi and perform the 
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FIG. 6. (Color online) The dispersing LLs of a three di¬ 
mensional WSM. We have chosen the following parameters 
A jt z = 0.2, tx,i/iz = 20, and B = 2T. The Landau levels 
with N=0, 1, 2 and 3 are respectively represented with the 
colors red, black, blue and green. 


Peierls substitution k —> k — eA/h. However, we do 
not linearize in k z , for illustrating how a dispersing LL 
behaves in the entire first Brillouin zone —tt < k z c < ir 
along the c-axis. In the Landau gauge A = (0,Bx, 0), 
the Hamiltonian becomes 


Hwsm ~ hv±[-id x (Ti + ( k y - eBx/h)a 2 \ 

+ ~[2t z cos{k z c) - A]<t 3 , (61) 


where v± = 2t±^a/h. The dispersion relations for the 
Landau levels are given by 


£iv=o(fcy, k z ) = -2 t z cos (k z c) + A, (62) 

(2 t z cos (k z c) — A) 2 


ew^oiky, k z , a) = sgn(a) 
+2 N 


h*vV 1/2 


1 2 
l B 


(63) 


where a = ±1. The dispersive LLs for N=0, 1, 2 and 
3 are showed in Fig. [6] for a particular choice of tight 
binding parameters. The LLL crosses the zero energy at 
k z c = ± arccos(A/2t z ), and the two minima of the higher 
Landau levels also occur at the same locations in the Bril¬ 
louin zone. An explicit consideration of Zeeman coupling, 
only shifts these minima and the Fermi points of the LLL 
to new locations, as long as (A + A z(B)) < 2 t Zl where 
A z{B) is the Zeeman splitting. However, this does not 
introduce any qualitative change to our conclusions re¬ 
garding the magnetoconductivity in the quantum limit. 
Both at and away from the half-filling the system be¬ 
haves as a metal and we can always find a set of Fermi 
points. If (A + A z(B)) > 21~, the minima of the N ^ 0 
, a = 1 and N = 0 LLs occur at k z = 0 At the same 
time the maxima of TV ^ 0 , a = — 1 LLs are also lo¬ 
cated at k z = 0. At half-filling the N = 0 LL is still 
partially filled with k F ,o — 7t/ 2, and the system acts as a 


metal. We also note that, for the N ^ 0 LLs, both up¬ 
per and lower components of the spinor wave function are 
nonzero and they correspond to and <pn- i, where p>N 
is the wave function for N -th Landau level of an electron 
gas. In contrast, for the LLL the lower component of the 
spinor wavefunction vanishes and the upper component 
is ip 0 . 

When 'J2Kv_l/Ib > (2 t z + A), both N = 1 LLs corre¬ 
sponding to a = ±1 remain well separated from the LLL. 
This is the simplest situation for analyzing the LMC. For 
this relativistic model the Fermi wavenumber is given by 

k F ,o{ B ) = J+2tt 2 Z!<5, (64) 

2c 

where 6 measures the carrier density as the deviation 
from the half-filling. In particular when, 


0 <l%< 


47tc<5 


< 


2 (hvF) 2 
(A + 2t z ) 2 ’ 


(65) 


only the lowest Landau level remains partially occupied. 
In this situation, the calculation of LMC is identical to 
the one for the layered metal presented in the previous 
section. The LMC due to the ionic scattering still varies 
as ~ B 2 . In contrast, the LMC due to the short range 
scattering is given by 


°*{B) = 


4 e 2 t 2 c 2 

hnfU 2 


x exp 


cos 2 (27r 2 1%cS) ( 1 + 


(l + 47 rl%cSy 


( 66 ) 


At half-filling, the zero range point impurities lead to 
a constant LMC, whereas the neutral Gaussian impuri¬ 
ties give rise to a B-linear LMC (when < 2a 2 ). For 
a system away from the half-filling, an increasing mag¬ 
netic field strength, asymptotically brings kp,o closer to 
the half-filled value 7t/2c. Therefore, for a finite 6, the 
LMC due to the zero range point impurities grows as 
cos 2 (Bq/(2B)) towards saturation. In contrast, the fi¬ 
nite range Gaussian impurities still lead to a B linear 
LMC (as the argument of the exponential saturates to 
7r 2 a 2 /c 2 ). In the presence of both ionic and short range 
impurity scattering, the LMR can initially decrease as 
B~ 2 before crossing over to the short range impurity 
dominated behavior. Therefore, for a pure WSM (with 
no carriers coming from additional parabolic bands) the 
LMR in the quantum limit can show upturn to become 
positive only due to an insulating ground state caused 
by the electronic interaction. But, the existence of addi¬ 
tional parabolic, bands will make the LMR to follow the 
generic behavior illustrated in Fig. Q] even without any 
broken symmetry phase. For (2 t z — A) < y/2 Hvj_/Ib, the 
LLs with a = — 1 remain completely filled, and do not 
contribute to the conductivity. The condition for quan¬ 
tum limit is then determined by 


2 t z s\n(2ii 2 l 2 B c8) + A < 


\[2 hv± 
Ib 


(67) 
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and our conclusions remain unaltered. When (2t z — A) > 
y/2hv ± /l B , LLs with a = — 1 start to become partially 
hlled. Such a situation can be addressed by using Eq. (fl9l) 
and a numerical evaluation of the transport lifetimes 

Ttr,N • 

Before moving on to the discussion for other physical 
systems, we point out the essential difference between our 
results and the ones found in the existing literature. To 
the best of our knowledge, the magnetic field dependence 
of the velocity vf, o and the transport lifetime Tb^(B) 
have been ignored in previous work s 18 i 19 i 21 i 35 i 46 . In these 
past studies the velocity vf ,o has been taken as the Fermi 
velocity of the Weyl excitations in the absence of a mag¬ 
netic field, and 7>,.o has been treated as a phenomenologi¬ 
cal constant (by ignoring the matrix elements of the LLL 
wavefunction). From this one obtains a B linear magne¬ 
toconductivity (following Eq. 1171) . Through our analysis 
we have showed that both vf,q(B) and Tb t o(B) have non¬ 
trivial magnetic field dependence, which in the quantum 
limit can not be ignored on any physical ground. 

For weak magnetic field strengths, the axial anomaly 
induced negative LMR of a WSM can also be understood 
from the semiclassical calculations of Son and SpivakJ^ 
This approach is valid in the absence of Landau quanti¬ 
zation, i.e., ui c t « 1. In this case, the axial anomaly 
of Weyl fermion appears in the semiclassical equations of 
motion (in the plane wave or Bloch band basis) through 
Berry curvature. In contrast, our calculations are per¬ 
formed in the strong magnetic field regime lu c t » 1 
(particularly in the quantum limit) in the spirit of origi¬ 
nal proposal of Nielsen and Ninomiya. These two theories 
can as such be considered complimentary since they are 
developed for different magnetic field regimes. 

However, there is a common feature in both theories. 
Son and Spivak do not explicitly compute the trans¬ 
port lifetime (in contrast to our work where we explicitly 
calculate the transport scattering times in the system). 
Rather they made an important and tacit assumption, 
that the intervalley scattering time (between Weyl points 
of opposite chirality) is considerably larger than the in¬ 
travalley scattering time and the latter determines the 
transport lifetime. In the absence of a magnetic field 
intravalley scattering completely dominate the resistiv¬ 
ity (i.e. intervalley scattering is unimportant in the sys¬ 
tem and can be ignored). This can only happen for long 
range impurity potential, since short range point impu¬ 
rities give equal intra and inter valley scattering rates. 
This requirement on the nature of the relaxation mech¬ 
anism in the work of Son and Spivak is precisely what 
we have demonstrated from a Landau level based cal¬ 
culation, without making any tacit assumption regard¬ 
ing the transport lifetime. We explicitly show in our 
work that long-range impurity potential generically gives 
rise to negative magnetoresistance (which is what Son 
and Spivak find in their semiclassical treatment) whereas 
short-range impurity scattering (not considered by Son 
and Spivak) by contrast gives rise to positive magnetore¬ 
sistance. Thus, it seems that the finding of Son and Spi¬ 


vak is a special restricted case of our general theory where 
we find that the magnetoconductivity is positive (nega¬ 
tive) for long-range (short-range) disorder potential. 

An important question is whether one can smoothly 
interpolate between these two (i.e. Landau-quantized 
strong-field and semiclassical weak-field) limits. The 
semiclassical methods cannot capture the physics of the 
strong field limit (as soon as quantum oscillations come 
into play), as it is oblivious to the Landau level forma¬ 
tion. In general it is quite challenging, but perhaps not 
impossible, to obtain such an interpolation starting from 
the Landau level basis. However, the interpolation be¬ 
tween the weak-field semiclassical regime and the strong- 
held Landau quantized regime is outside the scope of the 
Boltzmann transport theory used in both our and Son- 
Spivak work. 

If we consider TaAs as a typical example of WSM, 
magnetotransport experiments^ suggest that the quan¬ 
tum limit is achieved around B ~ 8T, which is reasonable 
for a carrier density n ~ 2.65 x 10 1 ' cm -3 . Therefore, the 
magnetic held range 5 T < B < 9 T for which a negative 
LMR is observed in this material in Refs. [74] and [75] does 
not belong to the semiclassical regime considered by Son 
and Spivak. Nevertheless, the LMC in this regime shows 
B 2 dependence, and a subsequent upturn towards posi¬ 
tive value. From the comparison of the quantum lifetime 
(extracted from Dingle temperature) and the measured 
transport lifetime, it is found that the transport lifetime 
is almost 100 times larger than the quantum lifetime^ 
This essentially points towards the dominance of the for¬ 
ward scattering mechanism due to the long range impu¬ 
rities. Therefore, the observations in this material are 
consistent with our theoretical predictions. 

We mention, however, a direct experimental way of 
verifying our prediction of an impurity-potential-range 
induced crossover between positive and negative LMR in 
WSM systems. This can be easily done by systematically 
introducing short-range (or long-range) disorder in the 
materials either through radiation damage or by actually 
incorporating atomic level impurities (or charged impu¬ 
rities) in the system, and measuring transport properties 
of different samples using different levels of controlled 
impurity incorporation. Such an experimental technique 
was instrumental in settling the importance of long-range 
versus short-range disorder in graphene transport prop¬ 
erties^. We urge similar experiments in WSM systems 
to investigate the importance of long- versus short-range 
disorder in producing positive versus negative magneto¬ 
conductance respectively as predicted in our theory. 

Recently a very large negative LMR has also been ob¬ 
served in the WSM phase of TaP between 5 T and 14 T 
(see Ref. [77b- In this range of magnetic field strengths, 
the LMC behaves as B 2 , which is precisely consistent 
with our theory (with predominant long-range ionic dis¬ 
order scattering). In contrast to the other WSM TaAs 
where the LMR shows saturation around 9T and a sub¬ 
sequent upturn, no such saturation has been observed in 
TaP up to 14 T. Despite possessing similar band struc- 
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tures, these two materials do have considerably different 
saturation field strengths (B* for TaP is yet to be de¬ 
termined from measurements at B > 14 T). This further 
vindicates our result that the saturation or crossover field 
strength B* is nonuniversal. As discussed above, our the¬ 
ory predicts that the introduction of controlled amount of 
short-range disorder in TaP should suppress the crossover 
field for the change from negative to positive LMR. 


VIII. DIRAC SEMIMETAL: TWO CONES 

When two Kramers degenerate conduction and va¬ 
lence bands linearly touch at isolated points in the Bril- 
louin zone we obtain a DSM. The quasiparticles are 
then described by the four component massless Dirac 
equations. The DSM phase with two Dirac cones^ lo¬ 
cated at (0, 0,±fco), has been experimentally realized in 
CdsAs^—, Na 3 B n n i 12 The pertinent Hamiltonian is 
built out of two copies of Weyl Hamiltonian H 2 in the 
following way 

Hd, 2 = H W sm ® t 3 . (68) 

For this reason, the quasiparticles’ velocities along the 
c axis at two cones have opposite signs. If we linearize 
the dispersion around the Dirac points we can write the 
following low energy Hamiltonian 

Hd ,2 ~ hv z k z a z ® 73 + Hv±(k x a x + k y a y ) ® 73 , (69) 

where ay’s are anticommutinng 4x4 Dirac matrices. The 
two Dirac cones have opposite chiral charge, which corre¬ 
sponds to the product of the velocities along three direc¬ 
tions. For a moment let us assume that the Fermi energy 
is pinned at the Dirac points. Then a relativistic field 
theory based calculation of the axial anomaly shows that 

5,1 + 5,2 = 0, (70) 

e 2 

5,1 - 2 = —Jr E ' (71) 

1/1 

where 1 and 2 denote the two flavors of the Dirac 
fermions. Therefore, the net axial anomaly (flavor sin¬ 
glet) cancels due to the presence of the doublet. In such 
a situation one may naively think that the axial anomaly 
does not contribute to the transport. However, a LL 
based calculation will clearly demonstrate that such a 
conclusion is invalid for determining the LMC. 

The LLs with N ^ 0 are now two fold degenerate, 
with the dispersion relation still given by Eq. (1631) . A 
new feature is the presence of two N = 0 LLs with 

ejv =0 (k y ,k z ,a) = sgn(a)[-2f~ cos(k z c) + A]. (72) 

Notice that the two N = 0 LLs levels have opposite cur¬ 
vatures. The dispersing LLs for a DSM with two cones 
are showed in Fig. [7] This type of LL structure after 
accounting for a shift caused by the Zeeman splitting be¬ 
comes consistent with the recent LL spectroscopy mea¬ 
surements performed on Cd 3 As 2 <^ 8 The magnetoconduc¬ 
tivity in the quantum limit is determined by the electron 
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FIG. 7. (Color online)The dispersing LLs of a DSM with 
two Dirac cones. We have chosen the following parameters 
A/t z = 0.2, t±,i/t z = 20, and B = 2T. The LLS with N=0,l,2 
and 3 are respectively represented with the colors red, black, 
blue and green. In contrast to a WSM, the N ^ 0 LLs are two 
fold degenerate, and there are two N=0 levels. In the quantum 
limit as denoted by the dashed brown line, the carriers in the 
two levels are respectively the electrons and the holes. 


like carriers from the N = 0, a = 1 level and the hole 
like carriers from the N = 0,ct = — 1 level. For this rea¬ 
son each level follows a formula similar to the ones for a 
WSM, where we need to use the appropriate electron and 
hole like carrier densities. We also emphasize that there 
is no impurity induced matrix element between these lev¬ 
els. For a net carrier density S measured with respect to 
the half-filling, the Fermi wavenumbers of these LLs are 


, 1 ( A 

fa, < , = -arccos( A cos (j ,p rf) 


+ a n 2 l 2 B 6, (73) 


and an increasing magnetic field tends to bring the 
wavenumbers close to the corresponding values at half- 
filling. After substituting these k Bja in Eq. (l48l) we can 
obtain the explicit expressions for the LMC due to each 
LL. We find that the LMC due to the axial anomaly 
is finite even if the Fermi energy is at the Dirac points 
(as assumed for a relativistic vacuum). The finite range 
Gaussian scattering gives rise to a H-linear positive LMC 
and zero range scattering only gives a constant LMC. If 
the long-range ionic scattering is the only source of re¬ 
laxation, the LMC will again grow as B 2 . 

These observations remain true when Zeeman spin 
splitting can be ignored. For a sufficiently strong mag¬ 
netic field, the Zeeman splitting modifies the dispersion 
of the N = 0 levels to 


e N =o(k y , fc z , a) = sgn(a)[—2t z cos (k z c) + A + A Z (B)\. 

(74) 

If, 2 t z < (A + A z(B)), there is a gap between the two 
LLs and the half-filled system behaves as an insulator. 
For finding k B ,on we need to make the modification A —> 
A + A z(B). Therefore, the Zeeman splitting has the 
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tendency to reduce the size of the Fermi wavenumbers, 
and the LMC. The asymptotic form for the suppression 
of the LMC due to the Zeeman splitting behaves as 


(A + A z {B)f - 

4 tl 


0(2i z — A — A Z (B)), 


where 0 is the Heaviside step function. When the trans¬ 
port is dominated by ionic scattering, we initially expect 
a positive LMC varying as B 2 , followed by its eventual 
upturn due to the Zeeman splitting effects. Therefore, 
the qualitative profile of the LMR will be similar to the 
one in Fig. [Q till A z becomes very large to cause an 
insulating behavior. 

The quantum limit in CdsAs 2 is achieved around Bq ~ 
43 T i 46 ! 79 This is in qualitative agreement with a rough 
estimate of the quantum limit using Ib < n 1 / 3 . If the 
density n s=s 1.9 x 10 24 m -3 is taken from the ARPES 
measurements^—, one obtains Bq ~ 10.3T. In contrast, 
if the density estimated from the the Hall coefficient is 
chosen, which can be a factor of 10 times larger than the 
previous one^, Bq is found in the range of 18T to 41T. 
At present there is no published LMR data in the quan¬ 
tum limit. However, a small negative LMR component 
has been observed on top of a very large positive back¬ 
ground even for much smaller field strength (B < 10 T)— 
Beginning around B ~ 2.5 T, the average value of the 
LMR shows a downturn due to the negative component. 
This quickly goes away if the angle between the applied 
current and the magnetic field deviates more than 4°. As 
a material, CdaAs 2 belongs to the family of gapless semi¬ 
conductors, and transport at low temperatures will be 
naturally dominated by ionic impurity scattering. There¬ 
fore, we predict that a systematic study of the LMR in 
CdsAs 2 and related materials up to the quantum limit 
will reveal a B 2 dependence of the positive LMC and its 
eventual downturn caused by the neutral impurities (see 
Fig. HD and the Zeeman spin splitting. 

Recently, Ref. He] has reported axial anomaly induced 
negative LMR in NaaBi, for which quantum limit is 
achieved around B ~ 4 — 6 T. The ratio between the 
intervalley scattering lifetime and the conventional trans¬ 
port lifetime for this material has been estimated to be 
around 40 — 60, which suggests the underlying relaxation 
mechanism due to the long range impurities. It has also 
been reported that the LMC varies as B 2 in a wide range 
of magnetic field strengths IT < B < 35 T. An increasing 
positive LMR background has been ascribed to a small 
misalignment between the current and the magnetic field. 
The results in the weak field limit can be understood 
within a semiclassical theory— But, in the broad range 
of magnetic field strengths 4 T < B < 35 T (quantum 
limit) the semiclassical methods become invalid. The 
experimental results clearly corroborate our main find¬ 
ing that in the quantum limit long range ionic scattering 
gives rise to LMC oc B 2 . 


IX. DIRAC SEMIMETALS: ODD NUMBER OF 
CONES 


A DSM with an odd number of band touching points or 
flavors does appear at the transition between the strong 
Z 2 topological and trivial insulators— Various exam¬ 
ples of such materials have already been presented in the 
Sec. |U A representative massless Dirac Hamiltonian for 
a tetragonal crystal with single cone is described by 

Hd, 1 = 2tj^ > i[sin(fc x a)ai + sin(fc !/ a)a 2 ] + 2 t z sin(fc 2 c)a 3 

+2f'[3 — cos [k x a) — cos (k y a) — cos(fc 2 c)]/3, (75) 

where otj = <jj <S> 13 and /? = <to < 8 > t\ are four an- 
ticommuting Dirac matrices in the chiral representa¬ 
tion. The Hamiltonian anticommutes with /3"/5, where 
75 = &0 ® 73 . The single Dirac cone is located at the 
T point k = (0,0,0). The Wilson mass term oc t' is es¬ 
sential to prohibit additional Dirac cones at other high 
symmetry points located at (7r, 7 t, 7r), (77 7t, 0), (7r, 0 , 7 t), 
(0 ,7 t, 7 r), (7 t, 0, 0), (0,7 r, 0), (0,0,7r). In the continuum 
limit this term causes a k dependent hybridization be¬ 
tween the right and the left handed Weyl fermions ac¬ 
cording to 


3 

Hd ,i ~ + /3 

j =i 


'h 2 kj 

2m± 


h 2 k 2 

H- £ 

2 m z 


(76) 


where v± = 2ta/h, v z = 2 tc/h, m± = h 2 /(2t'a 2 ), 
ro 2 = H 2 /{2t'c 2 ) . In the presence of the momentum 
dependent Dirac mass term, [ 774 , 75 ] ^ 0, and the con¬ 
tinuous chiral symmetry with respect to 75 is absent. 
This is reduced to a discrete Z 2 chiral symmetry captured 
through the spectral symmetry criterion {.f/ 4 ,/ 375 } = 0 . 
For the calculations of the Landau levels, if we ignore the 
higher gradient term such as /3 (kx — eAi) 2 /2mi, the 
dispersion relations are given by 


e N (k y ,k z ,a) = sgn(a) 


t 2 sin 2 (fc 2 c) + M 2 {k z c) 


+2 n 



(77) 


where M{k z c ) = 2t'[l — cos(fc 2 c)]. In contrast to the 
N = 0 level, all the higher Landau levels are two fold 
degenerate. The LLs described by Eq. m are illustrated 
in Fig. [ 8 ] 

An explicit consideration of (3 (kx — eAx) 2 /2mx leads 
to a gap between the two N = 0 levels, according to 


ejv —0 (ky , k z , n) — 


sgn(a) 


t 2 sin 2 {k z c) 


M(k z c) 


hui, 


c,-L 


1/2 


, (78) 


where w Cj x = eB/m±. The size of the gap is Hui Cy ±, and 
strictly speaking the half-filled DSM with odd number 
of cones is an insulator. The ratio of this gap and the 
bandwidth is roughly estimated to be 10 _5 H t'/t, and 
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FIG. 8. (Color online)The dispersing LLs of a DSM with 
a single Dirac cone as determined from Eq. G3. We have 
chosen the following parameters t±,i/t z = 20, t'/t = 0.4, and 
B = 2T. The LLs with N=0, 1, 2 and 3 are respectively 
represented with the colors red, black, blue and green. The 
Landau levels with N ^ 0 are two fold degenerate. In contrast 
N = 0 levels are nondegenerate. 


can be neglected for all practical purpose. The formula 
for N 0 are lengthy and not illuminating. For this 
reason they are not presented here. 

So far, we have not considered an explicit Zeeman cou¬ 
pling. Usually, the semiconducting materials have very 
large and anisotropic g tensors. In addition, the even 
parity and the odd parity bands (symmetric and anti¬ 
symmetric combination of chiral basis) can have different 
g factors. For a magnetic field along the z direction, the 
generic Zeeman coupling has a form 

Hz = A z,+ <2375 + A z,~ < 2375 / 3 , (79) 


solution of a quartic equation are tedious, and for this 
reason they are not reported here. 

By restricting ourselves to the quantum limit, and 
away from the half-filling, where one of the N = 0 lev¬ 
els (say a = + 1 ) is partially occupied, we can compute 
the LMC by following our general strategy. We use the 
following wavefunctions for the calculation of the matrix 
elements of impurity potential, 

•>Pa=+i =-^Hk z ),0,0,v(k z )\ip o , (81) 

^a=-i = [~v(k z ), 0 , 0 , u(k z )] (fio, (82) 


where ipo is the wavefunction of the LLL for a three di¬ 
mensional electron gas, and 


u 2 {k z ) + v 2 (k z ) = 1 , 
u 2 {k z ) - v 2 (k z ) = 


t sin fc- 


£N=o(k y ,k z ,a) - A z ,+\ 


(83) 

(84) 


These form factors modify the |{/ e //(2/c F) o(.B ))| 2 ob¬ 
tained previously from Eq. (E5l) and Eq. (051) by a multi¬ 
plicative factor 


i[u 2 (fc Fj0 ) - w 2 (A: Fi o)] 2 . (85) 

In addition we need to use kp 0 = 2ir 2 l^6, where <5 is 
the carrier density, defined with respect to half-filling. 
The short range Gaussian impurities can again cause a 
R-linear positive LMC. In contrast, the ionic impurities 
cause a positive LMC varying as ~ B 2 . 


A. Bii-^Sbj, 


where A Zj ± respectively capture the sum and the dif¬ 
ference between the Zeeman terms for the even and the 
odd parity bands. In the absence of the orbital coupling 
to the magnetic field, A z ,+ transforms the DSM into a 
WSM, whereas A z ,~ gives rise to a nodal ring. In the 
presence of the orbital coupling, the dispersion relations 
for N = 0 levels are given by 


£jv=o {k yi k Zl a) = sgn(a) 


hid, 


c,L 


t 2 sin 2 (k z c) + <^M(k z c) 
2 - 11/2 


+ A z.~ 


+ Az,+(80) 


The axial Zeeman coupling with Az,+< 2375 , shifts the 
N = 0 levels without opening a gap between them. In 
contrast, the antisymmetric part A z , -< 2375/3 tends to 
open a gap between the two N = 0 levels. Therefore, 
within the N = 0 subspace, the existence of a gap be¬ 
tween the <2 = ±1 levels is the most generic outcome. The 
Zeeman coupling also lifts the two fold spin degeneracy 
of the N ^ 0 Landau levels. However, the analytical ex¬ 
pressions for the N ^ 0 Landau levels obtained from the 


Recently, a negative LMR has been observed in 
Bii-zSbj,!^ We believe this observed LMR is consistent 
with our predicted generic axial anomaly in the quan¬ 
tum limit. In particular, the observed B 2 dependence of 
sigma(B) agrees with our predictions for the LMC in the 
quantum limit, when long-range ionic scattering is the 
dominant relaxation mechanism. Given that the mate¬ 
rial is a semimetal which is closely related to a narrow 
gap semiconductor, it is expected that transport will be 
primarily determined by ionic impurity scattering. In ad¬ 
dition our theory also provides a simple explanation for 
the upturn of the LMR beyond B ~ 4T as arising from 
the residual short-range defect scattering in the system 
as shown in Fig. [I] 

Bio. 96 Sbo .04 is a compensated semimetal, where the to¬ 
tal number of carriers in the conduction band is equal to 
that in the valence band. There are three Dirac cones at 
the inequivalent L points. The Fermi energy intersects 
the conduction band around these Dirac points, produc¬ 
ing three ellipsoidal electron pockets. These pockets are 
related by a rotation of 120° about the trigonal axis. 
In addition, there is an ellipsoidal hole pocket at the T 
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point, coming from the valence band, which can be mod¬ 
eled by a three dimensional parabolic electron gas. In the 
strong field limit, the positive LMR is mainly caused by 
the response of the hole pocket. The existing literature 
on magnetotransport and Nernst effect suggests that the 
quantum limit is achieved around B ~ 3T, when the B- 
ficld is directed along the trigonal axisSi. This happens 
due to the very low carrier density in the semimetal- 
lic Bii-^Sba,. The density of holes is estimated to be 
4 x 10 16 cm~ 3 . An equal number of electrons are dis¬ 
tributed among three electron pockets. A crude estimate 
using Ib < n -1 / 3 leads to the threshold for the quan¬ 
tum limit Bo ~ IT, which compares very well with the 
experimentally determined Bq. Due to the presence of 
anisotropic mass and g-tensors, the quantum limit for 
some of the Fermi pockets can occur at a smaller Bq 
when the magnetic field is applied along the bisectrix or 
the binary axes^. These facts indeed corroborate our 
claim that the experiments are mostly being performed 
in the quantum limit. The upturn of the LMR beyond 
4T is also in general agreement with our conclusions re¬ 
garding the asymptotic behavior of the LMR, which is 
illustrated in Fig. [I] 


B. ZrTes and HfTes 

In Ref. H, a positive LMC varying as B 2 has been 
identified for ZrTes up to B = 9T. For a very small held 
strength a positive LMR is observed, which has been 
ascribed to the weak antilocalization effects. In addi¬ 
tion, the authors have performed ARPES measurements, 
which suggest the existence of linearly dispersing bands 
over a significant range of energy. At present there is 
some uncertainty whether the underlying Dirac fermion 
is truly massless or possesses a tiny mass. Again, given 
the small size of the Fermi pockets and the low carrier 
density, it is conceivable that the quantum limit is at¬ 
tained for a weak enough magnetic Held. At least, the 
observed positive LMC varying as B 2 is consistent with 
our theory. Additional measurements are required to es¬ 
tablish the threshold value of the Held strength for the 
onset of the quantum limit. At least the previous mag¬ 
netotransport measurements^— suggest the existence 
of multiple Fermi pockets, and a carrier density for the 
hole pocket to nn ~ 6.5 x 10 23 m ~ 3 for ZrTes, which will 
lead to the onset of quantum limit around a few Teslas. 
If the field is applied along the b axis the threshold can 
be further decreased. 

The authors of Ref. [IB have assumed a constant relax¬ 
ation rate for the axial charge, and invoked the existence 
of a chiral magnetic current for explaining the observed 
B 2 behavior of the LMC. The argument is as follows: 
(i) the steady state axial charge arising due to the ABJ 
anomaly, causes an axial chemical potential /T 5 , which 
is proportional to E • Bp,, and (ii) there exists a chiral 
magnetic current 

jcME = H5B/(2n 2 h 2 ). (86) 


This formula for the chiral magnetic current has been de¬ 
rived in the context of the unbounded linear dispersion 
of the Dirac fermions in Ref. HB, and in this context it 
is a nondissipative, equilibrium current. For a bounded 
dispersion of solid state systems, there is no equilibrium 
chiral magnetic curren t 84 ' 85 . This can be seen in the 
following way. A straight forward application of the for¬ 
mulas in Ref. [HI leads to 


jCME = ^2{ip N \d kz H\il> 

n) 

N, k 


2 nl% 

*=> n. 



dk z deN(k z ) 
2i r dk z 


0(fi - e N (k z )). 


(87) 


For an unbounded linear dispersion one sets c —> 0, 
and the LLL then represents two disconnected worlds 
of left and right handed fermions with different occupa¬ 
tion numbers. The chiral magnetic current arises from 
these unbounded modes. For a bounded dispersion as 
we have shown for WSM and DSM, the left and the 
right handed worlds are connected through the bottom 
of the N=0 LLL (finite depth). Consequently, when we 
perform the integral over k z for a partially filled LL, 
we find the chiral magnetic current is proportional to 
ejv(&F,Ar) — £jv(— fci^Ar), which always vanishes even in 
the presence of a ^5 ^ 0. As far as the completely filled 
LLs are concerned, their contribution vanishes due to the 
periodicity of integrand. From a calculation of diffusive 
magnetotransport of a WSM in a weak magnetic field, 
a positive LMC proportional to B 2 has been found in 
Ref. 22, in agreement with the predictions of Ref. IH and 
Ref. |33. Since, the applied magnetic field strength up to 
9 T is not a weak field regime, further experiments at 
higher field is needed to clarify various aspects. 

Compared to ZrTe 5 , HfTes has a lower carrier density 
n/, ~ 6 x 10 22 m~ 387 ' 88 , which leads to a lower threshold 
for attaining quantum limit. Therefore, future measure¬ 
ments on HfTes should also reveal the axial anomaly in¬ 
duced LMR in the quantum limit, which should display 
the generic behavior illustrated in Fig. [T] 


X. CONCLUSIONS 

In this manuscript we have shown that the existence of 
a finite LMC or LMR for any generic three dimensional 
metal is purely a quantum mechanical effect and is a di¬ 
rect consequence of the one dimensional axial anomaly of 
the dispersive LLs (see Eq. m and Eq. (1191) ~). The oc¬ 
currence of the axial anomaly is usually associated with 
the relativistic field theory in odd spatial dimensions, in¬ 
volving massless Dirac or Weyl fermions with unbounded 
linear dispersion. Therefore, one of our main findings 
that the axial anomaly is a generic feature of an arbitrary 
metal in parallel electric and magnetic fields may seem 
surprising and counterintuitive. Generally, the Adler Bell 
Jackiw equation is derived by evaluating a triangle dia¬ 
gram, by employing a Lorentz invariant regularization 





20 


method. Such methods are well suited for the problems 
involving a fluctuating background gauge field, e.g., a 
non-abelian Yang Mills field, and the axial anomaly in 
that case is caused by the instanton density, which is 
proportional to Tr[E • B]. If we apply such methods of 
calculation for a system with a parabolic or a generic non- 
relativistic dispersion, we will not find the axial anomaly. 

However, the axial anomaly for the abelian background 
fields is somewhat subtle. For fluctuating electrodynamic 
fields, E and B are orthogonal and there is no axial 
anomaly. Only in the presence of externally imposed par¬ 
allel electric and magnetic fields, the axial anomaly be¬ 
comes finite. Since, a uniform background field strength 
is not a weak perturbation, many subtleties can often be 
missed out in a direct perturbative calculation performed 
in the plane wave basis. In the presence of a uniform ex¬ 
ternal magnetic field, any three dimensional metal with 
conserved electric charge demonstrates Landau quantiza¬ 
tion of the cyclotron motion, and the dispersing Landau 
levels provide the appropriate description of the vacuum, 
and also serve as the proper basis for perturbative calcu¬ 
lations. The Landau levels of a generic three dimensional 
system has one dimensional dispersion along B and an 
exact degeneracy factor eB/h. One dimensional conduc¬ 
tion channels in the presence of an external electric field 
do possess charge pumping between the Fermi points (at 
the intersection of the one dimensional dispersion and the 
Fermi level), and this is described by Eq. (JT]). In a one 
dimensional system, this is related to the uniform accel¬ 
eration of the center of mass in the presence of a uniform 
external electric field. This is the underlying reason for 
the emergence of the axial anomaly even for a three di¬ 
mensional metal with parabolic dispersion. When E and 
B are parallel, the net charge pumping is obtained after 
multiplying both sides of Eq. ©■ 

Quite strikingly, we have shown that the appearance 
of a negative LMR is not in any way tied to the existence 
of an underlying Dirac or Weyl band structure. Rather it 
is intimately related to the type of scattering mechanism 
present in the system. We have calculated the generic 
axial anomaly induced LMC and LMR in the quantum 
limit in the presence of short-range neutral impurities 
and long-range ionic impurities. We have shown that 
ionic impurity scattering gives rise to a large positive 
LMC (7 oc B 2 in the quantum limit (see Eq. IM1) , while 
neutral short range impurities and point defects give rise 
to a non-monotonic dependence on the magnetic field 
(see Eq. 12711 , which is dependent on the underlying band 
structure. However, in the quantum limit the combined 
effects of neutral and ionic impurities initially lead to a 
negative LMR which ultimately becomes positive after 
passing through a minimum as demonstrated in Fig. |T] 
and Eq. [M] 

We have also emphasized that due to the low den¬ 
sity of carriers in gapless semimetals, the quantum limit 
can be reached for a moderate magnetic field. As semi¬ 
conductors generically possess low carrier densities, with 
large amounts of charged impurities (due to n and p type 


dopants), they are prime candidates to observe the axial 
anomaly induced negative LMR. We have shown qualita¬ 
tive agreement between our predicted LMR in the quan¬ 
tum limit (see Fig. [L] and Eq. [Mil and the recent LMR 
measurements in TaAs, TaP, NaaBi, Bii_a,Sb x , ZrTe 5 
and a-(BEDT-TTF) 2 l 3 . The observed positive LMC 
in Bii-^Sba, and ZrTes varying as B 2 up to a nonuni- 
versal threshold is consistent with the LMC caused by 
long-range ionic impurity scattering. In particular, the 
quantum limit in NasBi occurs around B ~ 4T, and the 
positive LMC varying as B 2 is observed up to 35T. Cur¬ 
rently, this material provides the largest window for pos¬ 
itive LMC in the quantum limit, and it lends strong sup¬ 
port to our theoretical results regarding the importance 
of ionic scattering in gapless semiconductors. In addi¬ 
tion, we predict that the observed upturn of the LMR 
is arising due to the short-range neutral impurity scat¬ 
tering. We have shown that the quantum limit in HfTe 5 
occurs for a smaller magnetic field strength in comparison 
with ZrTes- Based on this we have predicted that the fu¬ 
ture LMR measurements on HfTes will show pronounced 
evidence of generic axial anomaly in the quantum limit 
and the LMR will display the generic behavior consistent 
with Fig. □ and Eq. [Mil . 

We have shown that in layered materials with high car¬ 
rier density, such that the underlying Fermi surface is a 
corrugated cylinder, one can obtain an effective situation 
like the quantum limit in a tilted magnetic field, when 
the tilt angle is in the vicinity of the Yamaji angle. Fi¬ 
nally, we have shown qualitative agreement between our 
predictions and the recently observed negative LMR in 
PdCo0 2 . 


Appendix A: Boltzmann equation for 
magnetoconductivity 


Here we consider the Boltzmann equation approach for 
evaluating the transport lifetimes. We restrict ourselves 
to conventional scalar potential disorder, which do not 
mix the orbital or spin indices of the eigenstates. The 
linearized Boltzmann equation for the transport lifetimes 
of different partially filled LLs is given by 


^2 W ( N ^ k y^ k z,a] N',k' y ,k' z ,a) 

N',k' y ,k' z 

Vl V'(fc') 


'l~tr,N,a( k z') 


~T~tr,N' ,a( k z ') 


= 1 , 


(Al) 


VN(k z ) 

where the impurity matrix element squared is given by 


W(N,k y ,k z ,a-,N',k' y ,k' z ,a) = ^-n i J ^|C7(q)| 2 

x | (N, k v , k z ,a\e^ r \N\ k' y , k' z ,a')\ 2 6(k y -k’ y + q v ) 
xS(k z - k! z + q z )6(e N {k z ) - e N >(k' z )). (A2) 


Here \N,k y ,k z , a) is the normalized wavefunction for the 
N the LL and it can be a spinor. In addition, a represents 
the orbital or the spin degeneracy of the Landau levels. 
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For concreteness we will consider the LLs of the three 
dimensional electron gas. In this case, 


I (N, ky, k z ,a\e^ r \N', k' y , k' z ,a')\ 2 = e~ u u' N ~ N \ 


xid~7:',Mi 2 


^min{ N, N } + 


min(Af,Ar ) v r (max{iV, N'} + 1) 


(A3) 


By setting N = N' we arrive at the general intra-LL scat¬ 
tering potential of Eq. (1471) . The potential in Eq. (l24l) is a 
special case, obtained by setting N = N' = 0 in Eq. m- 
For short range scattering due to neutral Gaussian impu¬ 
rities or zero range point impurities all the integrals can 
be performed analytically using 


where u = q F l B / 2 and L™(x) is the associated Laguerre 
polynomial. The Eq. Ell) can rewritten in a more con¬ 
venient form 


2nriil B 


E [ ^\Ueff( N , N '-,q z )f 

N' J 


Ttr,N{kz) 


v N'{k z + q z ) 
VN(k z ) 


Ttr,N'{kz + qz) 


1, (A4) 


where U e ff{N, N';q z ) are effective one dimensional scat¬ 
tering potentials obtained after integrating over q x and 
q y . The expressions for these potentials are given by 


\U ef f(N,N'- 1 q z )\ 2 = 


T (mm{N,N'} + l^ r £ 




2 / q±^B 


x\U(q±,q z )\ exp - 


r (max{iV, iV'} + 1) J (2tt) 2 Z| 
I j |A ? ' / —N\ f q± l 2 B \ |2 

I min(Ar,JV') 


B 


(A5) 


f 


ax ^ol 


\T a (T\\ 2 - r ( a + 2/3 + 1) (a - l) 2/j 
i py n r 2 (/3 + i) a a+2p+i 


X 2 F 1 -/?, -/?; -a - 2/3, 


a(a — 2) 


(A6) 


The integrals for Coulomb potential can also be per¬ 
formed analytically for small integers. Upon integrat¬ 
ing over q z we obtain the following M number of linear 
algebraic equations for the transport life times 


Ttr,N{\kF,N\) 


U e ff(N, N ; 2/eF,iv) 


| kF, 


N 


„ U e ff(N, N ; \kF,N ~ k F N > |) U e ff(N, N ; |fcp,iv + k F N > |) 


2|^F,iV' I 


+ E T tr,N' (\^F,N' I) 

N'^N 


N'^N i n'^N 

Ueff(N, N ; \k F ,N — k F N < |) U e ff(N, N ; |fcpjv + k F N > |) 


21 kF, 


N 


21 kF, 


N 


h 3 


2mriil 2 B 


,(A7) 


where M is the number of partially occupied LLs. If we 
just consider the LLL, the transport lifetime becomes 


7~tr, 0 


H 3 kF, o 

2mn l l 2 B |t/ e// (0,0,2fcF,o)| 2 


h 2 v F>0 

2ml B \U ef f{0, 0,2fcF,o)| 2 


(AS) 

I 


as has been discussed throughout the main text. The 
contributions of the short range scattering and the long 
range ionic scattering to the conductivity of a three di¬ 
mensional electron gas in the absence of a magnetic field 
are respectively given by 


a s (0) 

a c (0) 


.( 0 ) 


m 


mn s llQ 


A^Erf( 2 afcF( 0 )) - ak F ( 0)e- 4a2fc ^ (0) 


ne 2 K 2 h 3 k F (0) 
m 2irme 4 n c 


log 


4fc|(0) 


4fcf.(0) 

Qtf(®)- 4k F (0) + <?ff(0) 


(A9) 

(A10) 


which have been used in Fig. [TV] We have denoted the &f(0) and qrp( 0) respectively. 
Fermi and the Thomas Fermi wavevectors at B = 0 by 
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